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I. INTRODUCTION 


ee eee ee ee oe ee ee 


The purpose of the Advanced Tactical Inertial Guidance 
System(ATIGS) program is to demonstrate the feasibility of a 
low cost inertial system to be used in the Air Launched Low 
Volume Ramjet (ALVRJ) cruise missile for mid course 
guidance. Within the framework of this stated purpose lies 
the intent to furnish moderate accuracy ina strapdown 


inertial navigator with high reliability of operation. 


The strapdown inertial system requires a computer to 
provide inertial reference , hence the possibility of 
extending the computer's capability by installation of 
Kalman filtering aigorithms iS seen aS an area for 
investigation. Previous work (ref. 1,2) in this field 
indicates that the computational burden associated with the 
Kalman filter limits its usefulness when position updating 
Systems in the missile give highly accurate measurements of 
actual position. Most of the aforementioned computational 
burden resulted from the on-line gain generation required by 
a non-linear model within the Kalman filter. Hence if a 
linear model with sufficient performance were to pe 
incorporated and the Kalman gains generated off-line and 
Stored, then possibly the velocity estimation errors which 
are largely unaffected by the position updates could be 
reduced. : 


The purpose of this study was then threefold 


1) Test a linear model of missile dynamics for use as a 


Simulaticn tool. 





2) Determine the inertial navigator accuracy within the 
six degree of freedom simulation when a pure position reset 
device is installed which provides position updates at two 


points along the flight path. 


3) Determine improvements in missile performance if a 
Kalman filtering scheme were installed to estimate missile 


states between position updates. 


The means by which accomplishment of the désired 
purposes was obtained were various Montecarlo simulations 
utilizing existing data on the proposed inertial guidance 
systen. Extensive work,both in testing of physical 
equipment and in Simulation, had previously been 
accomplished by various departments of the Naval Weapons 
Center, China Lake,California. Hence accurate data as to 
component performance were available. These data were 
utilized to construct models of the components for computer 


Simulation. 


The simulation of a strapdown inertial guidance system 


requires the nonlinear computations relating observed 
accelerations to inertial frame coordinates. This is 
normally accomplished Within the guidance-navigator 


algorithm by a suitably chosen set of state variables aad 
their related non-linear dynamics. The essential idea 
behind the linearization technique used in this report is 
that the non-linear calculations relating accelerations and 
angular rates to velocity changes within the inertial frame 
could be accomplished upon observation of the said 
accelerations and rates and utilized as forcing runctions 
for a linear model of system dynamics. This corresponds to 
a free inertial system with observations physically aligned 
in the inertial frame of reference. Thus velocity changes 
in the inertial frame of the strapdown guidance system would 


be the non-linear combination of accelerations,angles and 





angle rates which are treated as inputs 


The model dynamics then are simple linear eguations for 


which Kalman filtering gains can be calculated and stored. 


The proof of the above linearization technique would be 
ln coOmpariscn of the existing empirical performance of the 
ATIGS system with the observed corresponding simulation of 
the system without Kalman filtering installed. Ref. (3) 
provides ample data of the drift of the ATIGS system as a 
function of time under actual flight conditions. fThe data 
are for a pod mounted version of ATIGS installed on an A-7 
seurcraft. Information from this report indicated that 
ground test drift of the system was on the order of 1 
Nautical mile (nm ) per hour under controlled temperature 
conditions. Under free flight test conditions without 
temperature control, performance was degraded to 4 nm per 
hour. The temperature instability was not incorporated into 
the simulation due to current effort to provide corrective 
measures within ATIGS. Hence verification of the model was 
assumed if Simulation indicated drifts of 1 to 2 nm per 


HOUL. 


Once verification of 1nertial-physical model was 
assumed, the next phase of observation of effect of pure 
position update was commenced. Ref.2 in an unclassified 
portion, contends that the optimal weighting of filtered 
position estimates and highly accurate measurement of 
position is such that the filtered estimates are ignored. 
This being the case, the computational burden imposed by a 
time varying Kalman filter may be unwarranted. The implicit 
assumption here is that the velocity errors incurred in an 


unfiltered system are not significantly decreased by the 


filter. Therefore the position reset feature would be 
sufficient to provide required accuracy at mid-course 
termination. This conclusion was tested by simulation of 





missile flight under conditions of increased noise levels 
within the ATIGS system and comparison with "normal" 
performance obtained,which allowed visualization of the 
Magnitude of end point error incurred under conditions of 


pure position reset and different noise level sensors. 


The final phase of study was the filtering of the sensor 
outputs to provide more accurate estimates of velocity 
throughout the flight. No attempt was made to filter the 
position update measurements due to their reported accuracy 
(sigma=50 ft). Thus any gains in performance would have to 
be from the filtered estimates of position after position 
update and continuous filtered estimates of velocity. The 
primary indicator of accuracy in line with ref.2 was taken 
to be cross range error and cross lange error covariance as 


a function of time. 
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The approach selected for isplementation of ATIGS within 
an actual cruise missile, was to employ a high-speed 
processor to handle transformation updating, earth rate 
torgquing and other minor tasks in order to save on time 
requirements on the more complex navigation and guidance 
computer. This peripheral processor would then supply the 
central processing unit(CPU) with the necessary information 
that it required to compute position within the inertial 


guidance frame. 


Thus the basic ATIGS unit involves ring laser gyros and 
accelerometers providing information to the peripheral 
processor wherein after suitable transformation, inertially 
referenced changes in state variables are supplied to the 
Main navigation-guidance computer. The main computer (CPU) 
then has an auxiliary input from an external position 
measuring device, the Microwave Area Correlation 
fixer (MICRAD). The MICRAD system is designed to provide 
highly accurate position measurements at two or three 
preselected checkpoints along the route of flight. These 
position updates would then be utilized within the CPU to 


reset the inertial guidance estimate of position. 


At present the only filtering system instalied is an 
application to the initial alignment scheme wherein a two 
stage initialization process is used to align the missile 
inertial frame with the parent aircraft inertial frame. 
Filtering is not used presently during midcourse guidance 
due to the position checkpoint feature and the short time of 
miaght. 


11 





A block diagram indicating proposed ATIGS utilization 
within the ALVRJ is shown in fig.1. 
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ii. UNERTDAL SYSTE 


Ae SYSTEM CONSIDERATIONS 


The inertial navigation system incorporated in ATIGS 
consists of the Honeywell GG-1300 Ring Laser gyro (RLG) and 
the Sunstrand Q-flex accelerometer acting as sensors. 
Empirical performance data for these sensors are found in 
table (1). Both sensors are considered to be of the 
integrating type in that the output of the RLG is in the 
form of total angle change per pulse and the output of the 
accelerometers are total velocity change per pulse. The 
measurement iS accomplished within the RLG by means of a 
counter system which totals the number of fringe pattern 
passages during each pulse period and similarly within tne 
accelerometers, the total velocity change is proportional to 


the magnitude of the output pulse. 


In view of the above characteristics it was felt that 
the inertial system as diagrammed in fig.4 could be modeled 
Simply and linearly by using the outputs of the sensors as 
forcing functions vice part of the state vectors. The 
Similarity tetween fig.1 and fig.4 should be noted. fhis 
would result in a net reduction in number of state variables 
by allowing the missile dynamics to consist of second order 
equations of displacement and first order equations for 


angular motion. 


Be. INERTIAL FRAME OF REFERENCE 
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1. 


i 


LVRJ Implementation 


The proposed ATIGS application to the ALVRJ utilises 
a local-level co-ordinate frame for navigation to the 
target. The local-level frame 1s characterized by 
North,East and up as the respective axis of calculations. 
The non-spherical nature of the earth introduces an angle 
calculation which relates the local gravity vector to the 
position vector of the origin from the earth’s center. In 
the ATIGS unit, the gravity vector calculation 1s 


accomplished by an inverse sguare gravitation nodel 


G =-( KM/R) R (1) 
where 
G gravity vector 
K earth's gravitation constant 
M nass of the earth 
R position vector from earth center to vehicle 


which approximates the local gravity vector to the desired 


degree of accuracy. 


The vector output of an orthogonal set of 
accelerometers is the geometric sum of all forces which act 
upon the vehicle and of course gravity is included. Since 
the above calculation is dependent on position, then the 
gravity vector is not constant during the time of flight. 
Thus to distinguish between the effect of external forces 
applied to the missile and the change in the gravity vector 


ah equation such as ees (2) 
eee Reo 
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a force exerted on instruments 


a coordinate transformation relating inertial 


axis(1) to accelerometer axis (a) 


R, inertially referenced acceleration 


G gravity vector 
resolves the time varying accelerometer outputs. 


ATIGS accomplishes the above procedure within the 
Missile and calculates the proper direction and range for a 


direct steer to the target. 
2. System Simulation 


The Simulation of the ATIGS mission began by 
approximating the local-level inertial frame defined in 1. 
above as a dcwn-range,cross-range and up frame of reference. 
Due to the limited range of the missile the gravity vector 
was considered constant and known, hence the simulation 
Simplified to a simple cartesian co-ordinate space wherein 
the navigator assumes knowledge of initial position, target 
position and range to target. The correct heading to the 
target was assumed to be the positive x-direction with the 


Cight-hand system defining positive cross-range accordingly. 


The initial position of the inertial frame of 
reference was taken to be the origin and instantaneous 
headings were taken to be the difference between the 


longitudenal body axis and the positive x direction. 


The ALVRJ maintains a constant wings level flight 


and this restriction was also placed upon the simulation. 
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C. GENERAL DEVELOPHENT 


Defining a set of state variables for an inertial systen 
undergoing arbitrary two dimensional translation and making 
the careful restriction that small angles and very small 
angular rateS are involved, one obtains (neglecting all 


noise inputs) 


if Distance from Origin down-range 

X Velocity in down-range direction 

se Distance from track centerline 

¥ Velocity component vertical to centerline 

Theta Angulaz displacement of body longitudinal axis to 
centerline 


Due to the small angles and negligible effect of angular 
rates one can approximate the accelerometer outputs as 
Z4=accel. in longitudenal body 


axis = Aip= AVid (3) 
fpeaccel, in Vaterat body axis= 
Aap= & Von 


where the B subScript indicates body axis. The output of 


the single gyro is 


2 = § *AT =486 (4) 


eee) 1 it 
Now 
a AV,,cos 8, - 4V,,sin 8, te Vip acos 8, 
“V5, Asin 84 (5) 
i= AV,,sin 8,+ AV,,cos Q, + Vind sin e, 


Woy Acos 8, 
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Where 


A X=Z,cos 8,-2,sin 8,tV,,4 cos 8,-V5,4 sin me 


4 Y=2,sin.0, +Z,cos 6, 7V,,A sin 04 V,,Acos 8, 


Which can be further approximated by 


AXE2, -2 -V,,8, & 8,-Vo,A0 


294 1D? 4 on 4 


$24 -2584 Voy 43 (7) 


A 22,0, 42, 1V,,00,-V5,8, 8, 


Be 457 4%3 


Further utilization of small angle approximation yields 


fabs X 

e cep 
Voy= a 
Hence 


AX=2Z. - _ 
2) > 259, - ¥2z 


AY = (8a) 
ea 25 + X Z. 


and for unit time intervals the discrete state equations are 
X(k+1)= X(k) + X(k) + .5*AX(k) 

K(kt1)= K(k) +A X(k) | 

Y(k+1)= Y(k) + Y¥(k) + .5* AY(k) (9) 

Y(k+1)= Y¥(k) +AYX(k) 

e{k+1)= Q(k) +46, (k) 
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Thus the observations can be treated as inputs to the 
system after appropriate substitution. The above model can 
be expanded to three dimensions and six degrees of freedon 
by the addition of one cartesian and two angular coordinates 
which would then consist of 


Z(k),Z(k), 51 


for a total of 9 states. 


This development was accomplished Without noise 
considerations. In the physical system noise would exist in 
the form c£ measurement noise in both the accelerometers and 
the gyros. Hence 


Zs=AV,, 4% 


Zo=AV5, + 8; (10) 
Z.-h0, + Y 

and 

a = a4 eC Lo Ye -¥, 9, -190 


? ° _ y* - 
2,9, +2,% x2. O,+ % +X QO 


> 
KX 
I 


The purpcse of this approach was to utilize the outputs 
of the sensors as forcing functions for the linear model. 
Hence in the preprocessor the non-linear calculations 
involving observations and states can easily be accomplished 
such that one then obtains 


? e 
U, = a4 — 259, = 127, 
t @ 
Un = 2,9, + 25 + XZ, GZ) 
? 
U, = Z 
E 3 
as hypothetical and known FOLCInG | func tons. Thus 
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Sewbstlitution into the model of the systen of U, 4) for Axtk) 
ema Ut) for Ay (x) would result in 


X(kt1) = X(k) + X(k) + .5*(AX(k)+ 4, 
+ v7.8, + ¥(k) PY ) 
(+1) = K(e) + AX(k) + 4 + 3, (k) 
ee) (13) 
¥(k+1) = ¥(k) + ¥(k) + .5*(AY(k)+ % 9, (k) 
+ ¥ +X(k)Q) 
Yktt) = (ce) +A Te) + 0, (ce) + ¥ 


+ X(k)Y 
6, (kti)= 6,(k) + Ae, (k) +¥ 


Hence the net result is the addition of a process noise 
term to the model. Analysis of this noise term proceeds 
With the systematic elimination of the non-linear tern 
involving x and © . This is easily justified due to 
the small value of & and © . Thus one is left with 
down-range process noise involving 2%, and YQ and 
cross-range noise terms involving oS and x 9 : Clearly 
the non-linear terms will dominate. The conclusion that 
logically follows is that the linearization technique 
utilized above will obviously be accurate for smali angles 
in a manner proportional to the magnitudes of the quantities 
YY and xp : Furthermore it indicates that Kalman 
filtering will be most effective in the estimation of angle 
information and of much Smaller benefit in the filtering of 


accelerometer noise. 
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PARAMETER UNIT PERFORMANCE UNITS 


RANDOM WALK (°/hr) GG-1300-RLG 


BIAS STABILITY (°/hr. ) 
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Figure 4 —- BLOCK DIAGRAM FOR SIMULATION OF ATIGS 
INSTALLATION 
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Thus the total state vector would consist of 9 states as 
opposed to the 15 state vector considered essential in 
reference 2. The above technigue was inspired by Kortum in 
his development (ref. 6) on Kalman filter applications. The 
inertial computation scheme is based on several assumptions, 
all of which are results of the short time of flight-short 
range requirements of the ALVRJ application. These 


assumptions are 


1 Constant gravity vector over a 40 nm. flight path 

ee Earth torquing not required 

ee Small angle assumptions for largest portions of 
maght 


D. BASIC SIT*ULATION DESIGN 


1. General Considerations 


In order to verify that the nine state system would 
be adequate for modeling purposes, a simulation program to 
test performance was conceived. The actual missile flight 


profile includes accelerations after launch 
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up to a maximum of 16 g's(fig.5)from the initial conditions 
given in fig.6. In addition, it is to be noted that for 
Simplicity a homing type guidance command is given to the 
Missile dynamics. This is recognised as ineffecient in 
practice but is simple in isplementation and provides a 
basis for comparison purposes. Since missile flight 
controls respond to airspeed and not inertial speed, the 
velocities used in the missile dynamics portion of the 
algorithm are true airspeeds. However the inertial systen 
must always compute in inertial velocities and hence must 
either adjust its calculations to include this difference or 
accept any error that this difference may entail. It is to 
be noted that no means of velocity measurement(l.e. 
doppler,mach gauge,etc.) 1s tec be provided. In this 
Simulation the difference is ignored due to the large 
magnitude of inertial velocities obtained and the short time 
Seeet light. 


2. Noise input Design 


The randcm number generators used in this simulation 
were of two types; Gaussian and uniforn. The Gaussian 
generators provided the noise inputs to the various sensors 
and the uniform generators provided the initial conditions 
and wind effects. The wind effects were such that constant 
bias directions of positive cross range and negative 
down-range conditions were imposed. The mean value of wind 
components in each direction was 30 ft per _ sec. with a 
range of + /s-8 ft per sec. maximum change per second. No 
attempt was made to ascertain the relevance of the chosen 
wind model,its purpose was purely to introduce a bias into 
the system equations in order that the scale factor noise 
term of the gyros could be exercised. The scale factor tern 
was finally dropped from the model of the gyros but the wind 


bias was retained. 
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The Gaussian generators provided noise inputs to 
each of the six installed sensors. The chosen model of the 


accelerometer noise term was 


= + *A (14) 

12,37 20G1 2,5 TH ,2,3%41,2,3 

where 

EQG random bias term held constant over the entire 


flight but varied prior to each sample in the montecarlo 
WG-scale factor term which varies through the flight 
The chosen model for the gyro noise term was more 


complex consisting of bias terms and a random walk tern. A 


random walk generatior is described in general terms as 


Eo=u (15) 
16 2 
where E. is the error at a given instant and U- 1s a 
white noise term with a standard deviation of GG. - The 
variance of E-U grows linearly with time according to the 
relation 
2 a 
Ur it Ou (16) 
r ie 


with Ge. given empirically in table 1 aS random walk in 
Yhr with an uncertainty of 191he =.0075 - For a sampie 


generator to be used every second 


G, = 5, / 60 7) 
a r 


An error from each gyro is introduced into the 


inertial computation which is treated as a change in © : 


() (t) = EO + G(t) (18) 


PHI error of each gyro per interval of time 


2/7 





EO constant bias term per flight 


G result of random walk 


3. Changes in Heading Resulting in Velocity Changes 


The computation of velocity in the inertial frame 
consists of terms which relate the change in heading to 
changes in inertial velocities. For small angular rates the 
changes in velocity in the inertial frames show as inputs 
1.e. AWXX,AWXY, AWYY,AWYZ,ANZZ, AWZX. 


AWXY 1S the change in velocity in the X(down-range) 
direction due to a change in direction Oy ~~ thas Ss 


expressed in small angle approximations as 





AWX 4 -A%, AB x 

y (19) 
AWY =| Ae, fl A8,, x6 
awz | Ae, - Ae, 1 7, 


AWX is the total change in velocity in the x direction due 
to small angle change. This calculation is performed in the 
appropriate missile dynamics portion of the simulation where 
the A@s are the results of commands from the guidance 
system of the inertial system. In the inertial system these 
quantities are treated as ASMS or measured Changes and all 


velocity changes are computed based on gyro outputs. 
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The specific algorithm for generation of guidance 
commands Was Simple due to the hcming type control employed. 
Inherent in the cross-range,down-range reference frame is 
the knowledge of distance remaining or “"time-to-go" for 
termination cf midcourse guidance and initiation of terminal 
guidance procedures. A parameter that continues to be 
Significant within this project is the small angle, small 
rate assumption. In the guidance algorithm the one_ second 
time intervals chosen for use would require an inordinately 
long sequencing operation if commands were given in terms of 
rates. To clarify this statement, rate systems require an 
initiation and termination command, which for: one second 
intervals would require a two second execution time. 
Therefore the guidance system employed within this project 
determines tctal angle change necessary and then commands an 
automatic pilot to accomplish this change. Thus the forcing 
function to the inertial navigator equations is not an input 
to the rate variables but rather to the angular displacement 
variables. No process noise was assumed for small angles 
hence the actual angle change was set equal to the commanded 
angle change. Notice that this is basically an open-loop 
process wherein the inertial navigator does not predict the 
next state based on the commanded heading change but rather 


on the noisy observed heading change. 


The logical guestion then arises as to the effect of 
system drift during guidance. Normal procedure would be for 
@a heading change command system in which an error signal 
generated by the navigator would be driven to null by the 
rotation of the vehicle. System drift during the heading 
change operation would result in process noise inputs to tne 
Navigator equations. The above was felt to be undesireable 
due to Kalman filter operation characteristics wherein 
Steady state gains are non-zero for a linear system under 
process noise. The method of circumventing this discrepancy 


was to use as the forcing function the observed heading 


Za 





change for the update of the navigator. Thus in the absence 
of process noise an accurate indicator of measurement noise 
would be the difference between observed heading change and 
commanded heading change. This concept was to be used 


during the Kalman filter application. 


The algorithm for command guidance iS given in 
fig.5. The "time-to-go" concept allows for a continual 
estimate of distance remaining, and a simple heading 
Calculation computes the heading change necessary for 
homing. 

I nel Gocition - present positicr 
sas position - present position (20) 
~ 20000-X(k) 

It should be noted that the fee profile 
Simulation was terminated at a point where the inertially 
computed down-range position was greated than/for equal to 
the final position. This would correspond to the switchover 


point for terminal guidance. 
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STATE VARIABLE COORDINATE STD, DEVIATION 





downrange 0.0 ft 
(0.0 m) 
POSITION crossrange O20 ft 


(0.0 m) 
altitude 35000 ft 
(10668 m) 


downrange 670 ft/sec 
(204.2 m/sec 


VELOC ITY crossrange 


670 ft/sec 
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vertical 0.0 ft /sec 
(0.0 m/sec) 
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Figure 7 - INITIAL CONDITIONS AT LAUNCH 
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5. Position Update System 
The position update system(MICRAD) is designed to 
fix the missile position at various check points along the 
meaght. Currently it is intended for the missile to 


havigate to each checkpoint 


its position. The missile 


the next check point and pr 


of 


CoO 


FOr purposes 


reduce complexity and 


missile estimated positio 


position at two discrete 
navigation 


SEC. 


The 
measurement system is @= 50 
Was added to t 
ef 


projected ac 


noise tern 
the 


empirical data. 


Simulation in an 


The magnitude of th 
1s the fundamental argumen 
the Kalman filter from the 
behind this 
are available (i.e. Kalman 
MICRAD then the 
would be in 


estimate,logically the 


assertion 15S 


estinate) 
f 
MI 


heavily 


two estimates would 


X=M X_ + (I-M) X 
Nn m 


then be 


inertially and upon arrival fix 
would then compute the course to 


oceed to navigate to that fpoint. 


this Simulation and in order to 
mputer time reguirements, the 
nm was set equal to its actual 

time steps. The inertial 


system was thus reset at time=15 sec and time=80 


curacy of the MICRAD position 
© ft at low altitudes. Thus 
he measurement of position within 
Lorre o with 


a 


retain agreement 


e deviation of the position fix 
t in ref.2 for the elimination of 
The 


that if two estimates of position 


inertial systen. rational 


filter position estimate and a 
weighting placed on each estimate 
avor of the 


CRAD fix. The optimal mix of the 


more accurate 


(2a) 
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fae Navigation estimate = X + p 
= navigation estimate error 


eemeasurement estimate = A +r 


m 
Yr = measurement error 
teem / (r + p) 

thus 


K =(r/(rtp))X, + (p/(r+p))X, (22) 


Since logically 
r<a¢ Pp 

X = X 
m 


and the Kalman filter estimates are ignored. 


From the above, it can be seen that filtering for 
positon is reguired only in the intervals between 
observations and that the optimal mix of filtered position 
and observed position reduces to the observed position for 


highly accurate measurements. 


6. tatistical Formulations 


= as eee ae ee ae ewe a ee ee 


The primary measures of systen performance for 
inertial navigation systems are mean of estimation error 
a. e. the mean of actual position minus inertially 
Calculated position) and variance of eStimation error. The 
mean of estimation error reflects the result of bias within 
the inertial system and the variance of eStimation error is 
an indication of system reliability. The simulation progran 


adopted in this study evaluated the mean of poSition and 
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velocity as well as their variance at each one second 


interval along the flight path in addition to the estimation 


error mean and vaLiance. The final position states were 
also computed and the nean and variance presented 
seperately. It was felt that this information could give a 


qualitative ccmparison of the missile performance with and 
Without the Kalman filter installed. 


The mean and variance eguations for each time 
interval was computed using standard Summation and averaging 


techniques. Thus 


o,f, > x; (23) 


with the inherent assumption that the relative frequency of 
occurence is analagous to the probability of occurence. The 


variance was computed similarly with 


s* = (1/(n-1)) 2 (x,-X)? (24) 


A 
74 
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A. INTRODUCTORY REMARKS 


Ref.12 discusses various aspects of the philosophy of 
Kalman filter applications in a very concise manner. Within 
this discussion the practical limitations of implementation 
are specified; the foremost limitation being that of a 
prereguisite knowledge of the exact Statistical description 
for each random signal within the system. This a priori 
information determines the degree of optimality of the 
filter. 


The filter is said to be optimum in the sense that it 
generates an unbiased, minimum variance estimate of the 
states of a linear system from some noisy measurement of a 
Subset of those states. The requirements imposed upon the 
designer are : exact knowledge of the system dynamics, 
covariances of initial conditions, and noise inputs. 
Departure from optimality arises when either estimates of 
the above guantities are used or approximations to the state 
eguations with lesser state variables comprise the model of 


system dynamics. 


Previous studies indicate that best results are obtained 


with pessimistic estimations of design parameters and 
linearization of non-linear systems if possible. The 
pessimistic estinate of design parameters reduce the 


sensitivity cf the design to deviations within the systen 


While linearization results in off-line gain calculations 
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which significantly reduce the computational reguirements. 


The linearization technigue analysis of this study 
indicates that the model of the system error terms are most 
critical in the estimation of theta. The cross-range and 
down-range position error being most dependent on these 
terms. Since digital filtering computational requirements 
increase roughly as the sguare of state variables, 
effeciency dictates that the number of filtered variables be 
Minimized. Therefore it was felt that the best usage of 
Kalman filter techniques would be in the estimation of 
theta. 


This evaluation is supported by various references (6, 8) 
most notably ref.6. 


Be. GENERAL THEORY 


Given a plant characterized by the linear discrete 


equations: 


X(k+1)= Det 1, k)Xk) + A (k+1,k)UCk)+ Wk) 


(2>)) 
ay = GC(k)X(k) + Vtk) 
where 
X(k) Column matrix of states 
D(k+,k) State transition matrix for time k to time k+1 
A (kr, Forcing transfer function . 
Ww (k) Precess noise term 
2(4) Matrix of observations 
V (i) Measurement noise 
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With the appropriate assumptions concerning zero mean 
noise terms and knowledge of the covariance of initial 
conditions, the optimal estimate of the state vectors at 
time k can te arrived at through suitable use of a Kalman 
filter. This estimate will be characterized by a minimun 
variance of estimation error as its criteria for optimality. 


The Kalman filter equations are: 


G(k) = Pk/k-1)C(k)° [OP O/k-1)C(e) + ROK) 7) 
P(k/k) = [1-G(k)C(k)]  P(k/k-1) 

P{kt+1/k) = @ (k+1,k)P(k/k) DB © (+1, k) + Q(k) (26) 
K(k/k) = X(k/k-1) + G(e) LZ(k) - C(k)X(k/k-1) 

X(k+1/k) = O (kt1,k)X(k/k) + A (k+1,k)U(k) 


where the notation (k+1/k) implies the estimate at time k+l 


given time k. 
mee oPF ECiFIC THEORY 


The current method of filtering inertial navigation 
systems is by building the filter specifically around a 
given error model of the gyro. The most commonly used error 
model 1s a Gauss-Markov drift modei described by both 


correlated ncise and bias terms. 


D =- (1/¥ )D+r (27) 


where D is the instantaneous value of error, VY 1s the 
correlation time and re is a bias term. Statistically 


the time function has been found to be roughly equivalent to 


a random walk model such that 
D=w+e (28) 
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where Wis a white noise term. Since the above drift is 
colored noise at has been handled previously by 


incorporating it in the state matrix such that it is part of 


the estimation process. The state equations would then 
become 

= % 

(29) 

D= wtr 

Z= 6,+ D 


Thus the Kalman filter would be designed for a 3 state 
vector vice the needed two states. Also the presence of 
process noise requires a non zero steady state value for the 
gains of all equations in the update portion of the Kalman 
iPeLter. This means that improper choice of parameters will 


bias the estimatior for all tine 


The overriding consideration of this study being the 
desire to maintain as few state variables as possible 


prompted a cioser look at the above method of analysis. 


The inertial navigator had been modeled as @ first order 


system earlier (sect. I) 


OI(k+1) = @I(k) + @OM(k) (30) 


however the Kalman filter requires an observation matrix and 
the above model uses the observation matrix AGM) as a 
forcing function therefore for filtering purposes the rate 


Variables had to be redefined and a new state added. 


OL, (k+1)= OL, (k) + OT, (k) + A OM 


(31) 
OL, (k+1)= OT, (k) 


The new state variable SI, 1s to represent a small 
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angular rate which should be zero under normal conditions. 
Notice that the forcing function is solely upon the angular 
position and not upon the rate variable. Now the 


observation is the output of the gyro which is 


AQM(k) = A@(k) + OL, (k)* At + (9 (k)* eee) 


where ~P is the noise term for the gyro given by 
p (k) = EO + Gk) (33) 
EO is a constant bias and g(k) is a white noise tern. The 


noise equation then for unit time intervals are 


AP(k)/ At = 20+ Chk) =AV(k) (34) 


and the observation equation is 


AOM(k) = AO(k) + OL, (k) + AY 
a) 
= O(k) + OI, (kK) + EO + G(k) (oe 
Now ASO is a known forcing function generated by the 


inertial navigator therefore a new observation can be 


defined as 


A@M(k) - AG(k) 
@1, (k) + EO + G(k) 


2” (k) 


(36) 


thus the discrete state equations are 


QI, (k+1) QI, (k) + OL, (k) + A0(k) 
el, (k+1) = ol, (k) 


C7) 
mek) = eI, (k) + EO + G(k) 


Now if the constant bias term EO is added to the ©1,CK) and 


redefined as eI, (k) _ eI, (k) + EO (38) 
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Then it can be seen that the Sa term takes on the 
Significance of an estimated drift and the state equations 


then assume the forn 


pe okt) 
OL, (k+1) 


1, (k) + eI,(k) + Ae(k) 
el, (k) (39) 


made) = @L, (k) + G(k) 


Where g(k) is white noise and the state equations are now in 
standard form for Kalman filtering. 


It was necessary to add one more state per gyro 
Simulated but this brings the total state vector to only 12 


which 1s still a net savings in state variables. 


A block diagram of the algorithm is given in fig.09. 
D. KALMAN FILTER RESULTS 


The generai Kalman equations of part B above were 
applied to the angular state variables in an off-line 
calculation by applying pessimistic eStimates of initial 


condition variance such that 


2 =03 
(40) 
2 =03 
oe = 1.000 x 10 
8, (0) 


and the variance on the meaSurement noise to be 


a ~ ) = 1.000 x io 


The resulting gains were found to be simple tine 
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functions such that at time T=KA t 


1-(2/k+1) 


G, (k) 


(41) 
1/k+l 


G, (k) 
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V. ANALYSIS OF RESULTS 


A. VERIFICATION OF RESULTS 


Due to the high degree of simplification of plant 
dynamics, the basic plant model is felt to be weak. The 
linear dynamics applied to the given g profile produced a 
trajectory similar to the expected flight profile of the 
ALVRJ but (due to its simplified nature) with several 
defects. The descent produced a net increase in forward 
velocity after level off which would not be the case for a 
true non-linear model. This defect can be accommodated by 
the inertial system and hence was not felt to be detrimental 


to the purposes of this study. 


Further work in this area would be recommended in order 
that basic defects ae the missile such as thrust 
Mis~alignment, process noise in control actuators and sensor 
misalignment be introduced into the simulation . A better 
tracking scheme could easily be instituted such that the 
guidance algorithn could institute a more efficient 


trajectory. 


It was felt that for purposes of this study the plant 
model provided a reasonable approach to ALVRJ simulation. 
The time of flight and velocity profile are within the sane 
general order of magnitude as the more complex simulations 
provided by ref.1 and correspond with physical intuition as 


to proper missile performance. 
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B. VERIFICATION OF ATIGS SIMULATION 


It was felt that the Simulation of the ATIGS would be 
accurate for purposes of this study if the error growth rate 
and the standard deviation growth rate incurred by the model 
were close to the observed guantities set forth in ref.1 and 
bert.3. Reference 3 stated that observed drift in the ATIGS 
test unit was approximately 1 nm /hr in ground test and 4 no 
/nroe in airborne test. The simulation results showa 
numerical average drift of 3.35 nm /hr. which was felt to 
be in the range of the actual system. Ref.1 indicates that 
the Simulation reported therein had a cross-range standard 
deviaticn growth rate of 1400 ft./min. between the first and 
second reset positions. The Simulation within this study 
had a cross-range standard deviation growth rate of 1459.4 
fo. ALN. 


From the above correlation in performance the study 
proceeded under the assumption that the simulated model of 
the ATIGS would provide a reasonable background for analysis 
of position reset and Kalman filter performance in an actual 


installation. 


Ce EFFECT OF POSITION RESET ON SIMULATED PERFORMANCE 


The unfiltered position reset feature of this Simulation 
had the expected result of drastically decreasing the 
variance at mid-course termination. The pure inertial 
navigator must contend with both initial condition errors 
and integrated boost phase errors which combine to produce 
large scale variance at mid-course termination. The initial 
position reset occurs after boost is complete and thus 


Virtually eliminates the initial condition and boost effect 
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On position. However, as expected, the velocity and angular 
errors of boost and initial condition still have effects on 
the final value of position at mid-course termination. The 
basic ATIGS navigator without position reset was found to 
have a radial uncertainty of 1608.6 ft. at mid-course 
termination (see table3). The addition of position reset to 
the basic model was found to reduce this uncertainty to 
466.9 ft. However these figures reflect the inherent 
accuracy of the inertial navigator with the very low error 
terms in the sensors. If the noise terms in the inertial 
navigator are allowed to have their standard deviations 
increased by a factor of 5, thus simulating a very noisy 
inertial navigator, the basic ATIGS model without position 
reset demonstrates an uncertainty of 8646.5 ft., and after 
the addition of pcsition reset,4128.0 ft. Thus it can be 
seen that,even with a position reset very close to 
mid-course termination, large uncertainty of position can 
accumulate due to the magnitude of the velocity errors which 


accrue throughout the flight. 


D. EFFECT OF ADDITION OF LINEAR SUBOPTIMAL KALHAN FILTER 


The Kalman fiiter proposed in this report had the effect 
of decreasing the variance of Simulation behavior along the 
flight path for ail tested situations (see table 2) . If 
cross-range standard deviations is taken as the criterion 
for performance quality, it can be seen that the filtered 
performance is readily superior at all noise levels. The 
‘“normal" noise level exhibits a radial uncertainty of 248.03 
ft. at mid-course termination and the X5 noise level 
exhibits a radial uncertainty of 1083.1 ft. The fundamental 
reason for this increase in accuracy iS shown in table 4, 
where the velocity errors at final update position are 


compared. The filtered velocity estimates, even at the 
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higher noise levels, are such that mid-course termination 


position is within much more reasonable bounds. 


The overall effect of the Kalman filter proposed in this 
study then is to reduce the end point variance of missile 
position at mid-course termination in a Significant manner. 
The unfiltered updates of position are seen to be valuable 
when used in conjunction with the Kalman filter but do not 
insure adequate missile performance if high noise levels are 


encountered throughout the flight. 


48 








“utw/3z €SO€ 
uTW/IF 6S°H99 
utw/3F OFVZOL 


uTw/IF 96°S66L 





utw/3IF SzZOl 
uTw/IF 7° 6SH7I 


HONVa-O HIMOYS 


YOUU eusts-T 











NOILVY TIVLSNI 
dO adAL HLIM AONVNAYOANAd NOTLVIAWIS-2 aTavs 


uTw/3F oggel | szu/uu cfeeyt J = xufm ogtz | 
ujw/IJ 9€ZE ay /uu /g° ry/wu 6€y" 


= = ~ EB OS = 








uyw/IF 986Z1 ay /wu 9°€8 


uTw/AF 8°6EvV ay /uu 7¢°T ay/wu 67g°T | 
umw/3F L7Ly1 ay/uu 6°/Ol Ay7wWU 79°89 
utw/IF 91Oyv ay /wu 67y°T | ay /wu €/°2 


HONVU-G HIMOUD HLMOUD HIMOYS 
wows ewsts-[ | YOMNT AONVU-9 YOWUA AONVA-Q 


zy/wu yE9°7B | 











"TVNYON 


GX 


" "LYNYON 


ASTON 





SOILY 


SOILV 


AdAL Nila 


49 





NOILVTTVLSNI 


dO TdXL HLIA AONVNAOINAd BNIOd TWNIA-€ ATAVE - Le eanbrtg 
IF 6°S8LS 53° L°0S6 13 6°SLOL IF C°LLI 













ava SSULy 















5 LE°6S IF L°O06T IF COLI IF OVSOVC 















IF3 1°89 IF G°BEC IF 61 53 G°ILT 








‘LVWAON avy $8714 


"TVNYON gear Tt 
Sia 


‘IVWYON SOILV 


ASTON AdAL Nie 














IF YEO VE IF LL OES IF O'L9 IF 000072 











3 G°LOIL IF 8°O6S1 45 L°T802 Te yOse 




















43 CL 69 IF O8L8EC 


a 


IF €6¢ 13 91L7e 



















13 7 CCC 5 8°97 IF L°LCS JF L°LOVY 

















VAG LLG 1 706 ae TF TOMOn Gt c 











IF HISS8L 755997 IF 6° EEI8 







33 L681 












IF L°CESl IF OLE6C IF 6E° ELV IF OC6LVC 












3 CLOSE 45 Y°SLY IF L°OEST 





IF 796 








IF 6°918 IF 7° 6009 43 CL LIC IF O880VC 










BUsT S-T 


sonvat SURG 





50 


ae 





type run ea neat velocity erro 
@ « 
< std. deviatio 


EL iGS NORMAL 110.4 ft/sec 
ATIGS X5 DO2,4/ #£6/ Sec 
Mi CRRD™ NORMAL 104.65 ft/sec 
AT EGS D& x 642 ft/sec 


| 


ATIGS , MICRAD 62.8 ft/sec 
& KALMAN 





| 


AT EGS AICRAD 262.47 ft/sec 


Figure 12 - TABLE 4-VELOCITY ERROR OF TYPE OF INSTALLATION 
AT FINAL CHECKPOINT 


au 





APPENDIX A 


RESULTS OF ATIGS SIMULATION WITHOUT FILTERING OR POSITION 
UPDATE 
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Figure 13 - SQRT OF DOWN RANGE VARIANCE CATIGS ONLY? 
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Figure 14 - SQRT OF CROSS RANGE VARIANCE CATIGS ONLY? 
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APPENDIX B 


RESULT OF ATIGS SIMULATION WITH POSITION RESET ONLY 
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Figure 17 - SQRT OF DOWN RANGE VARIANCE CATIGS WITH POSIT 
RESET? 
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Figure 20 - SQRT OF CROSS RANGE VELOCITY VARIANCE CATIGS 
WITH POSIT RESET? 





200 


0.58 2000£-03 


SIGT 


0.644 000E-03 


ATED 


62 





2 pe rs te 


eet Ne eer om 


VAR OF ERROR 


MEAN CF ERROR 


VAR OF TRACK 


MEAN OF TRACK 


TLME 


OMOM 
oo0oo 


OoOso0 
NMNoOWM 
mh Lyf 
nr he 
ODoOe 
Sanwa 
* ees © © © 


cCooe0o 


NIN TOO 
oooc 


ecooc 
-=Ore 
rood 
Donvc 
moOMmo 
DWMANMOO 
oo e« ¢ . 
Oooqooc > 
| t J 


OOO 
ooo 


aos 
WINN © 
=U 
Nes 
onoo wh 

THIMONOS 
ese 8 ee 


OCooo00o 


5550 93 


NM 
CoO°O 


Ooc 


roo 
anon 


1 
0 
4 
5000D 05 


Ore NOME 
°e* @ @ © © 6 


OoOOO0000 


ANE SIND 


HOE HC OE OC OC 


NNO 
oOoo0o0o 


Sooqgs> 
OoaVs oy 
we Oh = 
Orme a MOU 
TARanmorws 
WO NOI NID = 
®t &© @ © 
oOooo0°0co 
t tit 


OOONM™ 
oOoo00°o 


( 
oSoes 
omDoonnwn 
MO = TL 
UN ad PL 
ON OWN 
TaAanyo 
* © @ © © @ 


Ooo0oo0°o 


oO OK OC oe OC OC 


WOHON at 
Oooo0°o0o 


oooacoo 
ode LN SO OY 
MO MIA Toy 
ODormno— 
ONO O44 
FHIMNsS LT 
* ¢ ¢ @© @© @ 


OOD 000 


Paw 
OCOVvNU 


oOono0o0o 
MIWOO 
MONT © 
wor AMO 
Sh mun 
AN = HO 
see © @ ® 


oCoo000%9 


iy Om ety, ohm om 


ANON SIND 


HEE OC OE OK OC 


ONONNS 
OO0002D 


ooo vec 
CO) ond Pee LF) = OT) 
NOM OW 
Orns To 
=n 0m on 
TTMOag 
* ee © @ @ 


Cooooo°o 


OoocoD 


rrnaw 
WOUWO 


J 
onaoe 
owe 
WON tN O 
NNOG 
NSF OWN 
NAAN 
o¢ © @ & @ 


COOoo000 


ANIM FINO 


HO OE OO OK OE OC 


EE DS AP ee A ey ee a mh 


OTTO NINO 
oOOooo0o°o 


SS aoe 
Pe MD ey amt ON 
MHOnm 2 Om 
Ph Ph FOUN 
nt OT OS 
FAY 

*eeeees 


ooo00o°0oc 


NNMANO © 
Ooqoo0oo 


oDSa5 0 
NOM Or 
NaH Oa: 
TASKopP— 
O~OID Oe 
DW NIN OIIN OY 
eee ese? 
Coooo¢ 
t Pree 


0DOD0DOM 
Cooo°o0°o 


osoooco 
=“aADO@OWN 
WOO MN 
Fee ee 
ONODWN 
romano 
o © «© &© @ @ 


CO000F 


ANIM TI O 


WC DE HK OC OE OK 


Ol Om mMmry 
Coooo0o 


ooma0o&fh 
mOoOoorn 
m=ODTrOMmM 
PTA AH 
NO Omtaity 
TOM ot tty 
oe ¢ © @ ¢ 


Oooo°oo 


NAN OO 
oococe 


o50900 
PUD ONS 
ar iNaomr 
TOO 
DOr 
WON OYUN Dey 
eeeeee 
cooo°co 
I pious 


OO 0M 
oOoooo 


oooos 
HAM OW 
PLA OY AIUN 
SPINAL 
OO OWN 
TOSMNANS 
ove Foe 


ooo0c°oo°o 


Pr NAW 
OVvvuoo 


COoo.u® 
Onur © 
OaM~ © 
OD=—O00 
OAV nwN 
NN AMO 
es ¢« @ &© @» 6 
oooo°co 
tt 


te ey iy en glen 


aAalYO TINO 


FR OE HO OE OK DK 


OTOMMAN 
ooocoo 


Ooceas 
Nth rut 
CONN 
Nat a O OMT 
TW wasy 
TOAMaAAyM 
oe ee & © 


oOoo00°0cf 


NNO N= © 
OO000°0 


Soaqoegcs 
Wh Ooh 
MNO 
NAIM ~~ 
NINN OO} 
Pe OOO om 
see ®© @¢ © 
oOo0o000o 
! DL 


Om ONS 
oOooo0° 


Oo5050 
MN 
TM INDWN 
TOMMNNH 
OMOrTN 
FAIMING 
see @ @ 6 


oOo0oo0°o 


rIrNOW 
VwoOUUw 


OooooD 
ONAN 
NOD MO 

m= OONO 
NOANOW 
OnNSMO 
* © &® © @6© ¢& 


OO000CO 
t 


(ey Oe, ORS ot, Ae are 


ANIM TUN 


a? Se? te se te 


WE HK OC OE OK OK 


weir AOE Ger me ee SS ee OD et 


OTrOMMey 
ooO009%0 


oo°o0°o0o0o 
CWO Or 
Or HO 
OrmMmsr mm 
Pm AY CO 4 
PPro 
ee ¢ ¢« © @ 6 


ooco°ocen 


DO AOMMN 
ooo°ceo 


OSS 
NODaWN 
orm rn 
eat OY 
OWMWNOW 
TAM YO 
eee @ 8 @ 


oOoo00oo 


OTFOM EY 
OCcooce 


Sea eo 
OMOCoO. 
ATITOAWN 
Petr vod 
OM FING 
WF OY KY re a 

® © @ ¢ 8 & 


ooo00o°C® 


NAIM Aa OS 
Ooo g ce 


eoagacoe> 
aT ON OY 
Pe om FD OD 
OFTNOHO 
NINDS 
DONINNAIS 
2e® ¢ & 6 e e@ 
OCToo0o 
J Cp iat 


DO aAO MEY 
Cooo0oo 


ofD00°0 
Oman TWN 
MOoOOrWN 
MOS LN 
Cun nu 
TAMaAVo 
see &© 6 @ 


Oo0o000 


WTaOW 
VOOCOO 


oone eS 
JIMDMuUIh=w 
rTonmmoe 
Noaoe 
MIN OL 
anh Amo 
e © © © © 6 


oo0o00o0g 
t 


iy, oi, i, i, oi, 


ANIM TIVO 


ar Ge? se “er Ser to 


PC OE OE OC OC DK 


EE A tine ng tn ee mS 


63 





me ee ee 


VAR OF ERROR 


MEAN CF ERROR 


VAR OF TRACK 


MEAN CF TRACK 


Lobel's 


10 


Oz OM 
Qoocec 


CoOo00s 
AIAN O © 
ONIN OS Oe Ch me 
TIN oO 
LN FO) Be fre ON OY 
UN AT POON at oe 
eet @ @ @ 


oOoo00ce 


NNO Nae O 
ooocooen 


oOossoc 
Pe ot fh = ZO 
MOMmaONC 
INR Mem 
Mr NOS 
DONNNN ey 
eooeosve 
alotoje lola) 
| rete 


DO OM 
oooo°o 


C0o00a 
MAO ANN 
mFONIN 
Nn Orn 
Om st Cin 
Tas OANNS 
ee @ @ 8 @ 


oOoo0sco 


WS NOW 
wWwOUWO 


C5306 
WOC OW 
ING mM yO 
NNOOS 
5D Pew Pe PL 
ANIMA EO © 
ereeee 


COOo0COo°o 
I 


(SS Re ES Sp omy oe 
mI NM TST UNO 
eee SE eee ee te 


HC PK HE OK OK OK 


OTFTOMTMS 
oCCo°0co°o 


Oooo0oo 
NOOnwoOCt 
moron 
ON@DINOS 
Ose Rac 
DOF MANA wee 
ee 8 @ @ 8 


elojlelelele) 


NNN =O 
lekSola ele) 


aQereoeje' 
SBrAorown 
Om DAL aw 
ASN 
PY) TON so ome 
RRINN NN 
ee¢%o2¢@ @ 
2aO0000 
t itt 


Ct OM 
CO00oO 


Sesoac 
ce MON 
AN PR 
NW = Oin 
OnTnN 
FS me OOO 
eeest¢ @¢@ 


ooo°ocfte 


OFT OMIM 
Ooo0oceo 


O0OD00oO 
~-—O Om D> 
Qh P= 0) mt et WD 
OCD =O ON 
Ph FF Met YO 
OD POV OO) FO) ms 
oe¢¢ @¢ &@ @ 


oOoDco°co 


NNN =O 
Ooo0°0s0o 


200000 
DIAM ODN 
NON 
mIOOKS 
OrNrEmo 
ONNANNS 
* @¢e¢e¢ @ 6 
eo000co 
i tee 


Oo men 
COoOoCO 


S003c 
marm~on 
MONMIN 
PON we 1 
oOnm Or 
TF aMOANNO 
e*@e8 @ @ @ 


oOOoO00°o 


fee? See? See? Sge? tee Se 


OTrOMyrMs 
ooec[ec 


QOoong0eono 
aHMOONOM 
MAINE Nee 
orvoae- 
UIT COM COn 
per MM ON 
* @¢e 4 @ 8 


Ooo0o°eoe 


at NNNN GY 
eee eee 
oooeco 
t it 


OTOH TO 
CooCoo°c°o 


CQO00o0a 
fr Pm TP OA OD 
DO OP et hm ot wT 
NNO Oa 
MMO CA RAR 
aovrm®sms Ov 
eeeeee 


CoCcoco 


MANOA 
ols lojlelele) 


Gc 


Oe OMMS 
Ooooce 


ASQo00 
OTOM IN 
Ton 
TOV ON 
OUIN A an 
SAMS 
*>e¢©@¢ 08 8 @ 


Oooo°0eo 


WPM =i 
wOOOYU 


eaqaess 
OONOoO 
WNOUaO 
ONM—=0 
poe PLA ITT 
NN SMO 
ee ee ee 


OOOC0OVUSd 


Ay Cy in, iy ey om 


HIM TINO 


nent See? ge? See ge ee? 


OK OK OK OK OK 


Car INO ge er, 
ealeloe|s, 


QOooo0o0oo 
mNANOM 
ANOO +O 
—OCowomm 
Pe TCO IND 
TePK iN OO 
eo ¢6¢ «68 e686 


elejolejole) 


Ona VeaeO 
eleleleleo]«s, 


Heyjysoeyje' 
SDOmmMmN 
Uy met DL ome 
Wreas Oor 
OsrOov$d: 
ead NJ med FO) om 
e*e2ee8e ese f% @ 
OooCco°co°o 
1 | 


MO mim 
OWN OWN 
TOA VO 
es @e8 8 @ @ 


oCoo°c°esm 


WNT mL 
COwowo 


Oo0o0G 
wD0DAO 
YT DWIND © 
rN SO 
Om Oru 
NNSA OS 
eeeee @ 


ooo00°0 
' 


ine Gly PI ertiny oy 


MNO TINS 


wae? See? See? “a? nee ee 


PK HK CK PK OC 


TP MAM INE 
Oooo°o°0°o 


Qo0006 
INP OR TO OU 
TOM DU 
orornms 
Ware © N= 
Terni 
ee eet ee 


Ooo0e°ce 


om ON) eet ON) med od 


Ooceo0o 


I 
QOo0da04 
mont Ou 
met OCD MT) 
N~OrOoNn 
IVT NO iia 
INNO eh 
eee @ @ @ 
ooo0o°0° 
1 14 


OmOMey 
coo0e° 


Sg paeaeg 
Wer~oOown 
Pe ar ein 
mom Di 
OM Knir 
LPM OANS 
eof @ @ 6 


Oooo °coo 


WTO 
OVvoOOVU 


09000 
DONO 
TODOS 
LNOMO 
NT ore lta 
MANMOO 
ee «© @ @ @ @ 


COOC°0oOo 
| 


ANAT INO 


See? tee See See? See? “ee? 


PO PK HK PK OK OE 


16 


IN TMI 
elolejoajlele) 


Qo000o6 
Or Amn 
mi OO 
UNIO) md md oma CL 
Mm M—MOO 
amt oF ONL me re) 
o@e¢¢e8e¢8 


CoCce°oco 


lad hia © 
alovole lela) 


SoQqaqe 
HONS OD 
OAMNOM™ OD CO 
~NINAN ae 
MORO D 
NNN ao 
oe eee 
ooooc > 
I t t 


Om OMA 
ooo°0° 


COoa0adsd 
mt POLY 
MOm TI 
AD MO mith 
On— Din 
THIOFNO 
eee @ @ @ 


oooo9g° 


IAS MO mtr 
OVUVVUO 


OO3D60 
wWwvoUuImMoC 
TONAO 
=aNO0OC 
IN MPN 
MANNNMO 
e*¢ ee @ 


COoCoOo°oco 
I 


Ha NM FINO 


See? Seg? “ee? See Se? See 


Po KK OK OK OK 


17 


WT TMA 
CcCOoOOoCco°o 


a sogeas 
DDD ae ont 
SOmaomn 
OMNINNGS 
AMOWUOD 
TS SINS 
@ee@t@¢@e @ 


Oo0000 


NWNeANOO 
Cocece 


eOaagEeS 
SPO OM 
oma mMm&ry 
PO any Daw 
reno 
=NaIMN o 
*eeoee88°¢ 
ooo Ce > 
| tt 


Om ME 
oooco 


SSaae 
Torr 
DOIN DLT 
MNO®INW 
CM OW 
St ait erNS 
* ¢e@ 8 @ 8 


ooo°0cec 


LN mL 
VOCVVYYU 


OD200 
LrWUNDO 
SOONWND 
ONMrO 
mh OLN 
MNAYNIO 
e*eete 8 


Oooo0o0ce 


ee Os MEp oe a, OR 


HINO TINO 


eee Se? ee? Ne See? Se? 


PRK OH PK PK PK 


© 
_ 


INN = NOOO 
CoCco°oco 


ecaeaecss 
McOnae - © 
CNM mr ow 
OaNnoMNn 
Or Reno 
NN aN MS 
eee8e8e8 8 
oo°oc°oc°ocd 
I t 


O—OMM 
ooo°oc 


Q0n00 
Oman 
renin 
TOON 
ONO mIn 
PaM™TFNO 
e ¢@ e 8 e.|6°8 


oooo°0c°o 


OK OH OK OK OE 


oO 
—~ 


- 


= ey 


64 





rm re er 


a 


VAR OF ERROR 


MEAN CF ERROR 


VAR OF TRACK 


MEAN OF TRACK 


TIME 


OTL A 
SOcoce 


SqosSsyae 
UN =r ers OF 
SNOT O 
Awsome aw 
COWNDINOS 
=e OMA 
*eeere 


ooooce 


COMO 
ooeo°ceo 


OGgagage 
ATO rune 
OOS Os 
PD Pe ee UY 
Once nm 
TamNwow 
ee ¢ @ @ 8 


cosooc 


WS PO met TO) 
WwOUVOYU 


Seo So 6 
—~AOUMO DO 
STD ON 
NIA & uh 
nye omer TT 
SOI ene 
*eee?e 


Sec 


Fe mE ee, OR ED 
MNT UNO 
fae te? tee? tee tee 


IK OC OK OK OK OK 


oO 
A 


WTI NS 
cozece 


Sa coeck 
Nm~OaeNS 
Dam OOM 
morn 
WT OL eI 
ST tO SN 
eeeeve 


ooood0oo 


Oa ONO 
ooooc°o° 


Sooo qs 
One Om Cow 
Wey 2 OO 
rem Ou 
OMOWeL 
Pam IAN 
ee¢ @ @ @ 


oooo°ce 


WP Mm 
QDVwVvOoO 


OS000U0 


x OK OE OK OK OK 


_ 
™N 


OPW US 
eoqcec 


oooocs 
™OMNODOH 
Mr O eau 
PN OOF 
HSNO mrs 
NSN OUTEY 
eeeee?#es ss 


oo0oo0°ceo 


NIN BAIA 
oCceoceo 


Sao Ss 
WALA m AS 
rae oocm 
yeah CON oe 
DWM AOS 
MIS AYN 
ee @ee@e@ 
SceEaesS 
I It 


O~OMO Hs 
oooceoeo 


SOO Ja 
MNO 
Woo ..OO 
wen com 
OnOrm— wh 
THIMmNnNY 
eeveee 


ooo2too 


WT AS 


Ovovvow 


oOococs 
MO Wet WO 
MmDOOWE’? 
WN ONO eter) 
Mm Dame Ny 
TAOS Ot 
* ee @ @ 


a] he hem] er) 
1 ! 


ANDYTING 


>< OK OK OK OK OK 


ae 


WO FLAM US 
ele kale lle: 


oOnosSs000 
WON TIE 
MDMDMOOO 
Ore OWS 
COU UNA 
NIT MO-eN 
eooervee 


oDoooo 


NN ENN 
oooocese 


ogsse oe 
ONMODEe 
Co=aonrm 
Ur 0 Oop Cr 
NUNES Ot © 
TVS a 
e@eeee?eee 
Qooees 
' ted 


OOM Om 
Soeaae 


SOoDd0c 
NNDAM 
OOrms 
TOON OM 
OWMDAK~ WH 
SP mEPAL NAY 
* eee? «¢ 


opDaeooe 


WT OMS 
WOVOUL 


Sooo oe 
™ Orme OO 
Mom AUN 
MNS OM 
mr OM ON 
WIAs ey) = 
eeee# e@6hU6e 
cocecnd 

{ t 


le Oy ey OD 


SNOT NO 


wet tee Seer tee? ee 


re OC OOK OT OM 


m 
a) 


OP UIE OO 
oooece 


SConecoed 
Ooroor 
mah =r 
Mom ™MaW 
WUT NSS 
mrt att 
e® @ @ « @ @ 


Sa ceo 


NNN AIC 
cooo°0o°o°d 


Soca aS 
PANONGS 
PToo7RnnNo 
iON an 
Pm wh OFC ow 
PNAS 
eeee? » 
Sogo} 
{ (da 


DO mt OO ow 
oooceo 


SSooj\ ge 
WN ov 
Oomuny© 
Ton am 
CONnAKOeh 
STeMNY 
ee @ e@ 8 8 


ooocecot 


UVP OO mt WT 
COVUVOVUUY 


SScogse 
COODOO 
FO) OD met ome CO 
OoONoorm 
TARO 
WAY Tuy New 
eeeee#e@ 


oo°ooco°o 


iy Ry ey, OE omy, 
SAY UNO 
fae? tee? tee tee ee ee 


IE OK OK OK OE OK 


v 
A 


orwWmoy 
oooocesd 


Soocds Ss 
Mm Des ONY 
DVNNOrO 
AO OID 
TOU UV PA? 
TAT UT met 
oe ee re 


ooocooo 


NN aN OI 
o7UCcoceo 


eSeaeae 
Ors eo 
Sona 
Orman 
MAT OS 
WV OO PU ae 
oe eee 
eCeqoece 
I 1d 


Om OM Om 
COoocooco 


Ooosa5d0 
OmOm Hs 
AOoorm> 
PON DMomr 
ON ae ot 
TaATONN 
ee »« © @ ee @ 


oDoc°otzeo 


WT May 
OVVoUVO 


oooseo 
FTWDUINS 
MODwOwn Od 
mS OS 
Ooms DOr 
WII a Neat 
eeeee?e 


ooocoeo 


Pe ey OR OD OD 
mem east uO 
wae? tee Sar? ee? tee? ee 


I< OK OC OK OK HC 


uy 
tw 


OSI OS 
eee ace 


BOoOoOoOT 
OPKrMHOWST 
CnmOMNN Ss 
TORK OO 
MmoEoNOW 
WF OOD) 
ee6F ef @ 


Oooooe 


OM OL ow 
oooocoo 


(eleleyaja\(e} 
ros Ooms 
COoOmMr ND 
UNE Bh Oo, 
Cmastmm wr 
STaMonn 
eeee*?# e 


Sess So 


Wer gat 
OVvovvw 


Coo0enm 
MONTITO 
Mw rm 
YNONNM 
Om OIA A) 
Wes ON 
eevee ¢ @ 


ooocoe°o 


Te OD ED HD HD ay 


am srurd 


cae? tee? tee? wee ee? 


WK MoE OK OK OK 


0 
N 


OT UTE AS 
oqocece 


Oos9900 
ooanmmn”) 
MOM Os 


NN aN 
COOCCoSo 


eooaceqe 
OOODNAWH 
OOF MND 
Om ALO 
FIFOomnn 
WONT ON ot 
eee ¢ 8 (8 
Sete lSloye) 
I 1¢t 


Ot OO Se 
So qea 


Soosges 
Ono ar 2 
TINMOO 
MON OOM 
COWMNOMmLD 
PF eMoOona 
eeeetee 


MNejeis 


Ws (I~ WI 
OOovVvovuw 


DOQO000 
a OM OD 
MOUIMoOoo® 
aN TOM 
Nm Rw On 
Oosym Nee 
eo.°8 eee @ 


cCoooceo° 


ey ie OE OD ee, 
HAYM TF uVO 
eer ee tee nee te Se 


mK MK OOK OK OK 


~ 
NN 


ee a ES TE TE I RS is 5 A aD OE A 


ovromor 
oq oocae 


oOosgo000 
ammo 
STaNDOO 
samDooen 
NO OT uO 
mre qa 
e © @ @ ee 


oooec]e 


CWKOM Om 
oqCcecec 


eeyevere 
MAINONS 
LC ot = et CO 
NDDoom 
OWN OMA 
THUAN 
eeeete¢ee 


ooooeot 


WF MST 
OvVvovove 


esoagEe 
wororo 
NOOM Occ 
ONINE 
veer 
SO UT Qe 
eeeteee@ 
ocoocecde 

I | 


On om, om, Rom 
ANT UNO 
wae? ee wee tee? ee 


I OK OK OE OK OK 


@ 
A 


OoTFOMOS 
ovueeo°o 


SS5Seo06 
Wm COrn 
mm 
SAT COM math 
ha OD 
Ost a OMA 
oe ew ee 


oooooo 


ANIA NEON 
aoooccdwu 


Ssoegsae 
~eOnNone 
Oronmnw 
QnWDD> 
we MOA DON 
AST Om at 
eve e ee 
oooceo 
I itt 


WF AT 
WWOOVD 


Oooddo 
DM OmeMmro 
NDWNON® 
wy Nm 0 ue 
fe™~MOMnry 
WNNWNON a 
eeee?e 


ooo°o°eoe 
' I 


ie, Fin, Oy ay a 
MINT UO 
wae? San EP Seg? tae? “a 


mo OK OC OS OK OK 


oO 
N 


65 





Aah ee a 


ER a PE A OS 


Ce ty event te 


VAR CF ERROR 


VAR OF TRACK MEAN OF ERROR 


MEAN OF TRACK 


TIME 
30 


MmPromnor 
ooocoo 


SQooag 
Wh wr Greet NY 
DO met tN OH 
Orcr qr wv 
COOMNNr 
aT SOON 
eevee ee 


Coooo°ofD 


NN NM 
oooo00°o 


AQoo Qo¢ 
TlOAw 
MMO OO, 
Punt mene 
req OaeN 
BD NED) ot oe 
ee eve 8 
Sogo gre 
! et 


Det O 9 OD mt 
OoOcoooo 


ea_geaae 
COIrnowss 
Sat Da SO 
ODMH NM re, 
OWT Omi 
Path ONN 
see @ @ 8 


ococ]ce 


DPM ttn 
OCOOvOvod 


SaQDoao 
OON SHO 
NDNAW OD 
NWN DANE, 
oOn~mmm Ny 
Pee OUT GD ON ome 

ee @ @ 8 


e 
ooo°oecde 


ANNO rind 


aC HK OK OK OK 


her or ot 
Ocoocoo 


Q00000 
ot CN met ON LN ot 
SOC SHUM 
VSPNOLS HC 
=—300—r 
OS tt FO 
eee ee @ @ 


oceoc’e 


NN atanrmn 
Ccoooe°c°d 


SoockosS 


OmMONO 
ooocoo 


WAP AS 
VoOOOUwYYL 


eeaceeaS 
rouomM 
Na ~r emt UND 
ONO Om 
Ne Oey 
ANNO Ne 
oe © @ © ¢ 
Coovuceo 

t t 


ce iS, ite ee ee: 


SIN OAT UV 


ae? eee See? See? Se? we 


> OK oC OC 


= 
cy 


rroroy 
ooo°oce 


eqgocece 
PAWN WY ow 
NNN OW © 
QatrAaony 
MODOODD 
mt DP ment FN) 
eo @ @ 8 86 s 


ooceoc 


NN eNO AY 
ooo°o°o°o 


eseogage 
Pm Sun rin 
AM Snir oe 
NNO & & 
(AM a Omar 
COON et FT) eet et 
ss ¢@ 6 8 6 @hU8 s 
SO SOoce 
' treet 


DO mt OOS ot 
Cooeflz°o 


eSoHeaoae 
OUT OM 
INNO? O 
Ornom~mnm 
OWMOTm— ue 
See NM MUAY 
. ¢@ 8 s 8 s 


ofo000°o 


Nema 
VOUSOYU 


OOod020 
NOOmMUNO 
NO Nr = ® 
ONONGD FH 
AMO an 
P= WO Oh =e ow 
*e*eees 
oOoOoCooo 

! ! 


Sa NF INO 


KK OK 


A 
ry 


hmProror 
ooocceo 


O0Odb000 
at OT OO 
Hy DUN FO ome 
WOOD 
2SFPNOOWCH 
at PRN 
*® @ @©@ »« 6 8 


SoOoos © 


MN et NON 
Ooococo 


Oo0o000 
re ANNIS 
Or Oar 
DOIN OMDMNO 
Om sy at thy 
met ONY OS 0p eet we 
eeeee8 @ 
Soqoeqcc 
| ret 


Om ONO 
oOoooo°c 


eoaglEeoe 
AN NOW 
Tornoo 
Cm Ont a7 he ee, 
ONnM~Nr WN 
ees PT) OY TY 
ees eee 


MmOQOoo0oo0c 


BrP MAS 
VUWUOWO 


SQagagg> 
=O Tho 
CNC Ld hm P= CD 
MNODOWNM 
DM NOON 
Pee OY SO cet cant wd 
oo @ ee © @ @ 
ooupveco 

1 1 


me Re ee 
= NES UNO 
ce San? Sap? See? See? ee 


a OK OS OK OK OK 


mn 
aa) 


me ror or 
COoooce 


OOD00G 
LT me LT ap od 
rCOorome 
aA Mm ODO 
DOWN cry mt) 
ar Nat ON 
eee @ @ 8 


cooeceo 


Om OD oe 
oc oso © 


[an ia ea en = eat 
NM SPAS 
Ovrowoe 
OO ammo 
COonmornN 
Sf a Te NN 
sos ®@ & 8 8 s 


ooo lcs 


Mmyroror 
oo cece 


eagace 
RAN Hae S 
Weer com 
MO ONS 
Pe OD INI 
ae Nae 
oe ee eo 8 


ooocceo 


ANS MN 
co7p0°0e°d 


easaaqge 
COMNMNY 
ON AY 
Soorreys 
a SYN OR NI 
at NYS ot ot 
es es © @«@ © @ 
cococ]e 
\ tae 


Ome Oe 0 
e2oo00c°0o 


oo200C 
DOM Nm — 
rer ANNO CS 
noOo~raom 
OoCcorrn 
fare nn 
es @ @e@ @ 8 @ 


oDpoo°ce 


NF ANN YS 
OvOUHOw 


essaqaae 
za ons 
=—mononm 
rm OM am 
PP IAN ON 
DAO St at = 


elolelelole) 


Oe ome om OE 
MINIM INO 
a? See? ee? ee Se? 


Oe OS 


DO mt DOD met 
SCcooeoc 


Bya(sfa\eyje' 
OWNS Os 
NM 6 ONO 
Ose" 0M 
maONOR YD 
Sat TON 
ee e@ 68 © @ @ 


OoCCoCOD 


NASP MN AS 
OvVvooow 


COO0O00 
m-COMMNO 
=O hr =O D 
Frannm are 
ORM SN 
ODNOt at 
oe8e © @ @ 
COoOVvVOooo 

! 1 


m~ ow, SS Oe ee 
= NOFA 
cae See” eee ee? ee 


OK HK OK OK 


haw OT OW 
Socogcacs 


Qo0na00 
Mrmeiraer Oh 
DOD et oP ny 
POND OD et 
KHNANOCO 
AIS Oat ery 
ee & we © 8 


coocec 


mM rte ny mn 
Scoe sc > 


aoa eg> 
a7 an 
UN ON Pott ces CP 
KOS aN er) 
NIPNNON 
mt CD TO ot 
* es © « & @ 
oocoo°c]e 
' Vor ot 


Oe ONO— 
ooooe°coe 


Q0oca0agm0 
~morvor 
NEOs oO 
PON—Oom 
DFT O-MM 
Ta PT On 
oe ee & 8 


Oo tTOoD0eO 


WS MANNS 
Ovoovo 


ecassS 
OwOMwOUIO 
em OOS 
(UO OOM 
rh Denny 
ONY Ot et ot 
®* «© @ &@ @ 8 


ooo°oco 


wo, me 
AINA SrN Lo 


ae? egy? See? Se? ee? ee 


OK OK OK 


earwoes 
Sooo oc 


eoSseage 
FIONN D O 
MODOC s 
MD D=c 
MOM NOS 
OND eS eet et (T) 
ee e¢ @ @ &© 8 


oOoo°oo°c 


De ONO 
(Sole lersie; 


COD002 
r-ormnoy 
NO eS OO 
ODNODM™ 
AOI N 
Sie TENN 
ss ¢@¢@# &# @ @ 


Oco0°oD 


NAPSANAST 
WUMUYO 


eas sea qo 
TOQHD 
wt CO = CO OO 
DONNA AS 
a Own 
RNY Pe et et ot 
oeoe ee 8 
OoOD000eo 

! t 


a NO FNS 
See ot te ee ee 


> >t OE OK 


@ 
mr 


marOores 
cococcoc 


Ooo0000 
MAONO LN 
mr O™Os 
AM NINN CO 
WO mt ety at 
NPS ae 
ee @ es &© @ @ 


coooc°ce 


Ooaono— 
ooCceco 


Saga] 
NIN DONS ~F 
DOodWnMNMeo 
Dam OH OM OD Ar) 
me ODOM TA 
Pa TONY 

ees ¢ © @ @ 


oocoece 


AST ANAS 
WOOVOVU 


Sooaaee 
MO NHNO 
= sf Narco 
NANOS as 
Th WD ery 
ON ONY Pe et eet ot 
ee @#@@# @ # 8 
Se Ccocc 

| | 


SIAM S UND 


CK 


oO 
eC 


66 





———- 


VAR OF ERROR 


MEAN OF ERROR 


VAR OF TRACK 


MEAN OF TRACK 


TLME 
40 


we FU Ph 
oooo°oeo 


SQoaeoes 
CC) et et HP 
TOM NC 
=“MOOMCL 
CWMWNOMN © 
NSP St et et (O 
eeeee®# 


eoo000c 


MANNA AN 
ooeceoo 


Qo0o0oea 
DDOONC 
NTU et ™ ot OD 
OD et NN et 
MNDOMOLSMA 
ont CO) at OO} ON ot 
e@ ee @© @ @ 
ODoOO0O 
I ets 


OmMOMO— 
Qooocec 


Sou 
~reOnmas 
NNO OS 
omnd ond (FY TP 
—DDNM YS 
Tr me NCO CUT 
eooeevee 


oo°0°o00 


WYP MOAT 
OVOwWOO 


SSqaoqoc 
MON VUUNO 
—DAWNOnm 
NNT Qe es 
fmm Pm VU OP NY 
ONAN — Ore 
e@¢ 6 @ @ @ 


oo°T°00c 


at NOS 


OK OK OK OS 


mProerrs 
oooo°0? 


QGOD000 
TOmrmyns 
Tearoom 
STON WNS 
OOoomMmMean 
TAP at ome OO) 
@®eeseeeee 


oo000°0 


MNIANIN AN 
oOoo0ceoco 


eSoeeasgS] 
OOTmon 
Mm Oem Os 
mM Oorn 
Ii Qe OUTS 
emt NI ONQ OL ONT ee 
eevnvevee 
oo7Too°0 
| 4 G0 


Om OM Ca 
ooo°ec]e 


oaTouEe 
amt HIN OT 
MO meNNO 
ot ot HP Hiry 
a ON am LT 
wT et NO AISY 
eee2¢e¢*¢ 


oo0o0 D0 


NS ON So 
WOWOOvL 


Soacqgee 
CONNDO 
VCMmaOwD 


rroreys 
o00000 


ayelelelele 
TMDOMDO 
NO@aNes 
DOOR™OW 
NOT OSO 
CSU erty 
eeveevee 


ooooce 


OTMONS ST 
OVVIVHO 


OQoco0G 
at FO UU 
SO cet Pm ceed cee 
QO et OO) 
Om 0reN 
AN ema Re et 
@®eetee @¢# 


Qo0000 


mt NIT UNO 


wa? Oe ee Cee ee 


ah KE OK OS OC 


a 





~evyrores 
oocoeoec?e 


eSogaoce 
tate OU 
NOOUW TO 
TANKHOC 
WMO Tor 
PAO mt omy 

eeteeet se 


oococo°o°c 


ANAIN MAN 
ooo0oceo 


Sasaeasc 
TOOoNge 
SARANDON 
MN MIN 
mOOMmO™ 
mNN Te 
*#¢@ «© © # © 
oooo0oom 
I tte 


OT OMOWN 
cooeo°ee 


SSO 
GP HMO OAS 
ONO Or © 
aronDdy 
mat ot ONL 
STO UGN > 
eeee# oe ® 


ooco°oco 


OSMAN ST 
OVOOLO 


QO0000 
NWA LAO 


TAN SOD. 


WSN 
o~ mr AIC 
OT tO om 
* ee @¢@eee?é 


Qoo00e°o 


ANOS 


wa? eye Sa eee eee 


aS MK OE OK OK OS 


oe 





mTOT™ TS 
ololololeaia, 


eogace 
TAPIA WO 
onDynon~ 
meh COI 
COUUAS ™ 
TAIN OD matt CN] 
eee ¢@# @@ 


ooocec 


MNAIN MA 
Oo0co0°0o 


oastealde 
NNR IA 
Oa ™ 
SOODy Ss 
BONMOO 
emt md OY AY TNS 
ee @ @ @ @ 
coockzueo 
§ (yee 


OTM ST ST 
VUWOVULCY 


Coe sSeea 
mM ONO — 
CICA DO 
D at pe FF Ch SF 
Ondon~OO 
emt NS Pe et I 
* @ @ © @ @ 


Oo0000 
I | 


at NO PUY 
ea et COP Nee eee 


mS OR OK OS oe OS 


mr Owes 
oocec]e 


AVSlolSyleye) 
DOO AN amt MS 
Oovrronmnn 
oond Pom ont LO PF 
SPUN at ae 1 
eeoetvee 


oocece 


OPT MNSEM 
COOL0YL 


eZee 
,oorwrnm 
PaO 
O~ ON Mr 
OO OONT 
aad PN] Pom om 
*@¢@ef @ @ 


oOo0o000 


en Me Fe re a/~ 


mT TUN SO 


aq EO OS oe OS 


mPror™ Ss 
Doo°cDc 


ODOO000 
DOO uU, 
com mmen O 
art tS 
Pam FOU 
SF OP et et NY 
eet%eee¢e 


ooo°oce 


WNOMI™ O 
ODdecoeco 


eSeage 
NOMKOWNS 
LSNONMS 
MMOwOry 
Os o-rn— 
PD Dm a 

° eh. oe ¢@ e 


COoCO0Oo 


WPM M 
VVVO00 


SZaoogc 
WITTER CY) et 
CRD DD a 
PACD CO aN 
ont ot (OP LS 
amt YOO ee MIM 
ehmU6m©€¢ ¢* @e ¢ 
ooco0°o 

' I 


am NOS UO 


ee? cee ee eye Oe ge 


SE KK KOK 


~Toreys 
C$ooeo°0ce 


SSQoagge 
ond ond > me CH) od 
ar rN oe 
Tove ws Oo 
MAMAS ON 
VF ODMNN 
ee¢#e@@# 6 


ooeoce 


PANINI FO) ae 
SDOoocOooo 


Sac oas 
DOAN OW 
Am Some Oo 
NF OMNIS 
NIA S Me oO 
eeeeee 
OOo sag © 
I treed 


OU OM O 
ooo°oecoe 


eqeasasS 
GSDHDarTrwn 
NHoo™s 
toe OY OUT PD OY 
OMOA™~ a 
Ons OD ata 
e e 2 @ e e 


Seog. 


we? eee? ee ee ee 


mrorMy 
Ooo sD0°0 


Ses qqoe 
NNO SNe 
ODO SOOM 
DOT sv 
OOO UD 
WYP & at] a 
e¢#e# @#e 


co°ocece 


mPom™ 
ooco0oeo 


eooeed> 
Nga JF OI 
MOM MM 
WO OL 
OO =m OCONn 
at TM ONIN 
eee @@¢ @ 


Do°0°c°oo°o 


OSMAN TOS 
WOOOUO 


OO0N00 
MON NO 
&MAND orn 
CorOonm 
NAINA att 
m=) COON OT ot 
eee eo #@ 


oOoce0o 


mt IPT UNO 


>< OK OK OK OC 


mvroyses 
oocecko 


OQ00000 
ODOT 
FO me Mu 
Wt toon 
S$ MOMNO 
Wm QatQ\ an 
@eeee#*# 


cCoCcecoo 


Ph amt OM = 
ooocece¢e 


SeSeaqges 
Na aN UIT 
FMA FT OD 
Pe ry me UY UU 
Or meray 
amt SONNY 
eee # e¢e 


CcCooco0o°o 


OTKSMON SO 
OVOVOO 


SGOaooeds 
amt UN OP SP wnt 
OD D—DWO 
POY CN UN NO oe CD 
NIA ONS Om 
ot pr) ON MR 
eee#¢#e 


ocoo°o0c 
I 


SNM TF UNO 


mS MPC HK aK OK 


Ov 
+ 


2 ET EE Oe 


Oe AT 2 SN nat i IE tal na 


Gr 





VAR OF EPROR 


MEAN CF ERROR 


VAR OF TRACK 


MEAN OF TPACK 


TIME 


Pe ae Be Oe fh 12 


ooo°o°oec 


QOo00o0s 
DNOANS 
NNODITNH 
rTrOoOnoonD 
SnoOOoOW Oo 
ALY emt od J 
ee © @ @ 


ooo°coec 


MINION OO me 
ooo°oec 


SAaqoggs 
MAIO FS, 
SO a 
Soron=anh 
MDA oH 
MINIM EIN 
* © @ @ 68 68 
coo0oo0ce 
1 folie 


Pe od OP eet 
ooqooos 


aoc es 
TOMS O 
hm SNDNYY 
Oo~onon 
as tyh= Orr 
oe ¢ &© © © 


ecoooco°o 


OTN Oo 
WOUWUOOOo 


O2060 39 
rm FUN SS MOI 
DOMND WO 
OOD mi 
NN OS 
mt TCD EO) mt 
*eoe8 @¢ @ 


Oooono0oeo 
' 


ny ty, oy Ay ey ee 


Hae TUN 


EK OK OE KOK 


ere ore 
oooo°o0°c 


Ooooaao 
rrr ore 
MDH am eH 
THOmMmon 
CEO et ON 
WD Pe et et A 
oo e © o @ 


oOoo°o°oo 


ONO a 
ooo°o°0oo 


OOR008 
th Ay OD 
SNS MONS 
NOUS TM 
MINIM 
NNOMmMm mG 
*o ee &¢ 8 @ 
elololale ier 
I teins 


P~ O00) Pe 
Oooo°e°oo 


OOSD000 
NONSOM 
Net Ors cv 
Mm SDM 
Om Doorn 
atin mh CUT 
oe @ © 6 @ 


oOoo°eec 


ww TONS O 
COwowvd 


OO0Onaceo 
MANU D 
=“ON ODO 
ORO th =& 
TIALS ASD 
=D CUE) wet 
*¢@:¢ 8 © @ 


ooo°0°eo 


ie ing ay, Py ie 
SCY PUNO 
Sone? Ce eee “eye ee 


HR OK OK OK OK OE 


| 


he wh eh 
elelelolale: 


Qa0D00 
mor~onmnm 


MANNA MO 
COOoo°o°oeo 


oocsaoa 
OOrmnocsc 
Donnan” 
MONICIM a 
MM Pre wm L 
MNO MON 
ee @ @ & @ 
ovoose 
I tees 


frm wed oD) 10) he om 
elololelele) 


COGQO0D 
SOY eat CP eed Pome 0) 
DOM DO~ + LO 
f= O ah WY 
ODavon 
YOON 
oo © © @ @ 6 


oooo0°oeco 


OT MVS 
OUwooo 


QOB0000 
CUS O 
MOWED O 
MrT Of 
FA CIUYUAS OF 
=O A 
*ees © @ 


oOooo00e°o 


FT etm Sy ong oy 
HIN MFO 
wa eee ee ee ee ee 


HK HK OC OOK OK 


™N 
wa 


Pare ernm ys 
oooo00°ce 


OOd0o0o00 
COVUKnase 
MOC ONO 
Oaw>rasrn 
OMN OWNS 
Pee fee ound eed ON) on 

eo |¢ 8 ® s ® 


oOoo0000 


TOE CO 
ooooee9ee 


fre aD Ne ee 
oooc[e$9o 


OCcoooes 
nom wsaonm 
Wrm~Owr O 
we CO OO WD 
ODMmMenOrn 
=SIYODNAN 
oe ¢6@h6thlU 


oOoo0o°co 


OPN O 
OCOveVwosd 


Oonogng0co 
Wt OCR © 
WOUND O 
OOH Ohh =o 
MUU SF ON 
wt VOD OS PO) ome 
*oe © @ © © 


oOooe°o°oo 


ey i, a ay 


aN mT UO 


cane Wee Se ee ee ee 


HK EEE OOK 


aS 


Iw somone 


ooo0o°ce 


alelelele ye) 
=DONMOS 
MNO mM CY 
NM NOOO 
WMO OD 
Pee Pee cand oad J on 
eee @ 8 @ 


o°o°oo°oc 


POO ODEON 00, ome 
COOoo0vo 


Coosa 
O~=—OoO-— 
OPMNUIA 
SrnNnoac 
MOA MM a 
NA mmoy 
° s ® ® s s 
ooo°o°o°o 
; yo 


Phat OM 
ooo0ceo 


oOooaed 
OmsrAnm 
wONOMM LO 
mPONA SO 
COON O— 
emt CO DON ae 
oe et oo © 


oOoDoo°00o 


VPNs = 
OWwOuUwo 


QOooos 
KIWI VED 
POM Nno 
COD FO) et et DO 
MN OOS 
KINO 
ooeeeve 
eleleleleie) 

t | 


Fa a iD oat wate oR 


SANTINO 


ce? See? eee ee ee wee 


mS HK PK IK PE OK 


vt 
Ta 


Salar eas 


Ooooco0 


Sala lalel> 
TNVNDON 
OrnDvoamn 
ODOOITND 
aotrsr Oro 
Pee fe weed wend OJ ed 
oot @ @ @ 


QIOOOO0oO 


TAN NIC Miae 
oooc°oe°oc 


SoQqoIao¢s 
MOoOoCO ND 
OO a Neat 
UVP Oat 
MO OMOO 
NN Mme 
oeoeee¢ @ 
oooo°ee 
1 1iet 


MSF OM he 
ooooc°oec 


CO0000 
MANN =O 
he mt OP 
Mm OOM MO 
ON OOON 
=NVOUMNA 
* © «© @ @ @ 


oDoCo°oce 


we Moye O 
COOoOoouUY 


OCcoo0oo 
MmODDmmuwr 
mt YOO et 
MAM Oa 
TNyOOTO 
SAMOA 
[ee ° 6 


oo°o°0°e 
' 


il ey aay ay, a, 


SOM Tn O 


cee? eee” ee? ee ee? 


PS 9K PK OK OK OK 


2¥e4 


Pear he orm as 
coceco 


Oo3000 
Mon eyun DO 
MOOD Oy amt wt O 
ODO RON 
Teor oOooa 
CD weed ed OY ed 
® * ¢ ° @ ® 


oooo°c°o 


MINN NE 
ooooceo 


AaoOoaogaek 
MOTNOW 
Re mt OD OO 5 
Or ON 0 LO 
mmowraw 
NN mm 
® ® s ® tJ e 
Soong s 
1 red 


he SF DOMMme 
Qqoooece 


ODoDS000 
mm et met (OS 
(0) wet SE OIE CO 
OODORO SOD 
OND DO— 
SIND UN = 
*e* © © © © 


oooo00eo 


OP MOAI a 
OOWOowD 


Ond0oo00o0 
YP AO) 
TOOMAWN 
OKrOMODO 
wr or for ty 
SIMO AIM 
ee © &@ @ 


OOo000°0 
i ‘ 


AOI FUND 


wer er ee wer ee ee 


Hh KK OK OK OK 


Oo 
Ww 


Tm ome 
Oooocoo 


Ooo00na0 
OMY ONO 
heh er Or h 
TMD syy 
Hrawna Om 
COR tt mt 

oeoeeseoee?es®@ 


oOoo°o0°oo 


my om™iu 
oooceo 


Oooesae2 
OmNNAMS 
ANOS MS 
CORT I 
COr~nonon 
ATOMS 
o¢¢@ @ @ 8 


ooooce 


OTOP NY 
WOOOWSO 


Sa0000CS 
Omori 
~~oomm 
OO =a OOr 
SPO DMs 
SAY CO CU a 
oe eevee 
oooooeo 

' ' 


=m rund 


cae wee wee We ee 


Hh KOK OK OKO 


57 


Pe ame re 
oo°oo°eso 


Oooagoc 
™=OOMreS 
WNIND™ TH 
IN Sh CAS 
TWOOOS 
Om aan 
of © © @ © 


Oooceo 


rr omte 
SooeSse 


oooo°o0cdc 
ANAS OM 
O=amm om 
Oa S 
ConNnuoMm 
am O CUA ON 
oe. © @ @ 


oooo00°o 


OTMaNST O 
OOVOCVLO 


OD0000 
- tA OO 
WeerwOoneg 
NONNO HO 
Wea fre Hee om 
m4 10) CO OU OO) met 
eee? 6@ @ 
oQOvuce°oo 
4 


foe Deny any any ah oo 
aN rin O 


cee? Nee? eee? ae eee? ee 


PK PK IE KE PK OE 


58 


MNT mm 
ooooc]e 


oooaoacsS 
NTN 0 
TIM Nam OD 
am “Doc 
STOW aD 
NAO mn 
oe*eeoe¢ ¢ @ 
Oooeo0o0° 
I toe 


hm OT 
200000 


O235000 
UN he hm fh mm 
Maye am 
CS et NY et mt OD 
=sOrCan 
AAAI 
eee? @ 


ooo°0°0ct 


Or mr n 
VOvo0o00 


O20000 
US) od FY ot CD 
NPROnn 
0oo-mDrY-= 
Wom Om 
OOM 
ose @ @ @ 


ooooee:e¢ 


ele Clay Py ly ee 


=A TUN O 
wee wer wer er? er er 


PR PK OE OE OK OE 


Oo 
wn 


68 





em ese 


MEAN GF ERRDR 


VAR DF ERROR 


VAR OF TRACK 


MEAN QF TRACK 


TIME 


60 


Oreres 
ococoece 


Ooooo°o 
md LL) et FF) Pe ST 
NO mm Hr 
NM FEL 
ovrwonn 
wate mt PTT OY 
e e e@h6¢ e ¢ 


Cooooeco 


FLO) me ut, 
ooocoeo 


OS0902 
ONO 
CmDnon 
Mm UNO OO aw 
SLANT mt NI 
atee DINAN 
> ee @ e828 8 @ 


ooooce 


OPTMON SN 
VOOCVYY 


SosS oC 
Mymm™ WI 
WNT ONO 
KaON orm 
Noam NCES 
ot 1) OD PO IT) met 
eee? 8 @ 
oo°o°oeoco 

j t 


Pe tty ee ty oy 


HNO TUN 


jee See? See? See te ee 


HOC OE OS Om OE 


Qo me arm wr 
Cooe°co 


CO0O0S0 
ASoOMns 
eA mo 
Oe~mm~ wor 
=m KOM 
eat Pe et ONT ON 
ee ff &@ 8 


ooo°oeoce 


VOOs+O 
NN AOS 
e*e¢f8e 8 @ @ 
ooooeoce 
I roa 


LUND MO) Pe et 
o2mooeco 


Sooss4e: 
ar OnmM 
HSMM SH 
Tm ONIM ty 
AIMOwoNm 
wt mt CALA INI ON 
e*@e¢ 8 @ @ 


CoCoCoe°e 


Or MNT O 
COVOvLY 


-COOG000 


== Oh GIANT 
DWAIN Moo 
NNEOMrF 
OMNOMME 
mM yh ae 

e* ¢ e¢ @ @ e 
ooooeco 

{ 


at YEO TINO 


2S OK OK EOE OC 


61 


OrT~sres 
Oooeoce 


ooqgoe 
Qatar OU 
Wah OCR 
ODDO 
weDOOrm 
wnt me EY YON 
ee ee ee 


ooooc°oo 


MON MAN a 
OoooecSe 


Ooo0sdco 
TNHN Om 
DO CO awh hm O 
ONOMeVwory 
FT ODOmMOonm 
NA a OTOH 
ee ee ee ee 
oocoocesc¢ 
' ee 


miINO Mm in 
ooocec’e 


OOd005 
WT DON 
Mm SNES 
re NiO 
MPM Own 
amet CALA NINN 
e*¢e8@ © © @ 


oOooeooe°e 


VION 
OeVOVGWO 


OOoO2D000 
HOrIUCm 
at OOOM 
ONT SPE 
OWNOS OS 
mE EDM ww 
ee ¢@ ¢ © @ 
ODo000o 

| t 


ee Ss SE sy Oe 
KI NATO O 
wee? ge Se? ge Cee ee 


WS HC OE HE OH OE 


ez 


Ostresms 
COoooeeceo 


GOo0o000 
NNO 
TO HNSr 
NNAROts 
Ch wt eet OY 
ome BYE POEN 
eee @ 2 @ 


OoODoooe 


PIA A OO aw 
ooooes¢ 


Oo loeeac 
DOASNWNS 
WOONOWNOD 
mOor on 
Toon atrt 
NNSA re 
ee¢e? @ @ @ 
oooo°coe 
1 or on 


rmUNOM Mm IN 
oooeceo 


CO0O0000 
Age Ph LY et SD 
ml OOS aw 
WD OW at aw 
Toe OTN 
SON OIN NIN 
e¢ @ ¢ 6¢@ ° ¢ 


Cooeec¢ 


ormosr ty 
OVoooo 


Sosa saoe 
OmwMureryy 
Th Oem 
NO ume 
DO MInINE = 
mt PO) OD OT) et 
eee#8 * @¢ 
ooooco°o 

| 1 


a ety ae te oy 
HON STUN SO 


x OO OKO 


63 


Oty ohs 
oooeoce 


eoocege 
ot Pe LY he FO) 
TINHONM 
Om mt Srl 
NOME MS 


mr NIN 


COoooce 


MON MYA) a 
OoOoooc°oeo 


Sea qgae 
mA OU re 
TOO oo Ph fh at 
DaeOomsnh 
PO aS ath 
NN sre 
eeerv ee 
oooooeg 
! on 


mINOM~ UY 
ooocoeo 


Oo00sda0 
Om hw 
CONOR sO 
SRrOT Ss 
Om Dome 
aN OIA KIN 

2? 2 @ 


2 2 ® 
oocecso 


OTN STN 
VYOVOS 


OOO0OCOo 
ONMNONWOS 
mTOomOr 
NOND awh 
MmAUDOM ST 
mt) DD ae 
eee @ @ @ 


oooooce 


HN SUVO 


fee? tee? See? tee? Se ee 


XC HS HK HE OC OC 


64 


Owe om 
ODo0oco 


PNA me 
ocooce 


oo0cognsndo 


mUNomrm in 
oo°coeo°o 


ocecsee 
at OO DON 
Os ws OW 
WD Crm rm © 
Dew CO en ty 
at AMINAM 
eeo68f © 8 @ 


ooocece 


OTONSN 
CooVv°cVvd 


Sogoqcea > 
mii Oy 
EC at MO 
OTN ON 
mM MT DOU 
KHIM COMPO 
ee e@ © © @ @ 


ooeo°ceo 
! 


ee ey, Oe ee, oy 
NIT) PUNO 
fe? eet Se? te? eee ee 


MH OE OHO 


65 


cmd tnetath annithaies dontterlciahieetatiedtias 


wwe re 
ooooeo 


OOCAaOo 
CATO AU AUMm 
mOMMDN 
ONT weet 
ToNnmMos 
mere NON AT ON 
2 °° @¢ 2@ 2 ®@ 


Oooocoe 


MAN RY AU aw 
oooocoe 


Sooggc 
PO AY MDL a 
STNUMOMAN 
at OO M~FN 
LV OR at HO 
NN a Oe 


“ooooce 


Lets ae tat 


Pe LN OM in 
ooocec 


COoOoCcOoOOo 
DOMINO” 
ONIN =O 
NOWROUNS 
ae OTC 
NMR CAM 
eee eee 


Ooooceco 


ores 
OOoooeecoe 


C0O00O0nD 
NOOWM 6 
NDP MAO 
MmeAO OWN 
ST fh 0) te 
ate KIO ST ON 
. oe | ¢ 0? )68lhlU8 


ocoooo 


TOYO ON et 
oooeoo 


Saya 
NOM WO 


rm IN ORME UN 
Coooeo°o 


Cocooc]Nd 
D~ANOWN 
™ DMN OW 
Wat PT et 
TMM UIIIN GY 
NIT OHDOON 
ef @@ 6 @ 


opDooosd¢ 


OTN SEN 
OOovoow 


OC45030 
WrioO Murr m 
Oa NOmry 
NOSTNN 
DMMmyNes 
AOD SMO aw 
eeeree 


ooocoe°o 
1 t 


(A ny elt, OS ty ett, 


SNOT UNO 


we? Ge? See See Se Se 


OO OK 


Orne wes 
ooecoce 


eEeaqggaqge 
ooows,:re 
oOorwvrory 
mrs 
INH a wen O 
me ON ST ON 
e* 2 ® © @ ®@ 


ooooce 


TY NIE ANT) met 
COoooo°o 


Ooocoeoe 
NN at DP 
neo oP 
wt OO OD eet MOLT 
LA CR ONLIA RM ory 
NN et OS TF 
° 2 2 © @© ¢ 
oooocoo 
! trot 


MUA OM IN 
Ooooocoe 


oooocka 
DMM MNS Oo 
Owronms 
Dem DAIS 
DOOD 
NS DON 
ee ee 


Coogee] 


OFSMNST 
OvVUVovoe 


OoOosaooo 
Wma 
RON a UN 
UI Ome PT O 
DMMHeMAley 
A OOP Mh 
ee ee ee ee 
Ceooo0oce 

t | 


PR ett ety Ste ame, 
HAN TUDO 
Nee wae Mae Nee Ne tet 


HM OE OO 


68 


Were Ss 
coooeo 


SLsocqcec 
= DOM DM 
Oor~ nos 
OmMom De 
ODOM 
ate ONION AT ON 

* ° e e e e 


Oooooofm 


MAU Ne 
ooocec 


eecsaoee 
MPN SUN 
Pr NIU DOIN 
NOMM cor 
WV DT mS 
NNO re 
eee 2 @ 8 
oooooo 
t hota 


MmINOM ML 
ocoooeoc 


cooocca 
DPMMNOM 
WWM OR Se 
NSP IS 
(YD iw OOO 
MPT HM OT 
eeevree 


cooceooo 


OTM N 
OoOvuovecsd 


OOOO Od 
TP FT OC 
nmr © 
Arm sg tuey 
CD Oy OR UY ON mee 
attr oy 
ee 
Soooe > 

I | 


ae sty Ee tty fy, 
KOM TINO 


mh HOC OC OH OC 


69 


69 





1 erp mee 


EI a crea a 


_ 


VAR CF ERROR 


ERROR 


MEAN CF 


VAR OF TRACK 


MEAN OF TRACK 


TIME 


Ow>rmwzsess 
oooo°odc 


ooeoecde> 
NDDMNOD 
ON OO DO 
Peat ML FO 
O Nat Om 
tO MN FN 

eee? 8 


oooo0°0e 


CYNON EO met 
Ooscooos 


eCcsesgEe 
MmO Muro 
WOM at PO 
YF DONe 
Worm. 
Nye Me 
* @ @© © © @ 
ooooco° 
I { | al 


mA Mh 
oOo7mzOo7ToO 


eocsoae 
Tht Ounn 
momnurt Oo 
eit & OD 
Oa Ore 
em OO MAM 
eee 8 @ 


oocoocd 


OFM N 
OvVoovev 


oeosaeqge 
Mgr wre oO 
DWI URRY 
NDOAOr 
CMD OM 
mt TP MAY 
°° e6e@«f @ @ @ 
Ses 
' 


(ly i, Re gm, a 


at NON UNO 


fe? ee ee Se ee ee 


mM OC OK OK OK OX 


© 
ft 


Dare rr 
oOooo0co 


ooocoesS 
WONONTH 
DO=NAMO 
TO OUI 
MmONINO O 
at SOMONE ES 
e# @® $6 © @© @ 


O2oo0000 


TINO AN Mone 
OOO00°0 


AaoDSeaa 
Pt OMe 
AAS NO TD 
ITATHOOM 
WV Tee ey 
NVA OD 
e®e¢ ¢ 8 66 
Goo cec 
1 teat tie 


ROO MRI 
oOo00°0°m 


Soagee 
ODD mi Or 
anya mOny 
oor~7nor 
OD) at ont ] CH, 
PM OON S 
e#o © @ @ @ 


oQoooestoo 


OFS ON 
VOVUVWOO 


Sessa jjg 
PNA OY 
OD mt VIO 
OD atari fh 
SPO O CN ant one 
mt he FT wy 
oe0e@f @ &@ @ 
Cooooco 

| ‘ 


NOT UNS 
ee ge? ee ee ee ee 


2h MK KOK OOK 


md 
ft 


ores sr 
oo00°0c 


eoce sid 
TWH Mr 
NOW a 
NJ at UN OD om &D 
Oorroom 
“= DMS 
eee © @ @ 


ooos0co 


TAO OAS met 
Ovcooo 


Ooooce 
fre Ds et ome Pee Pe 


NAN DSO 


rWOTE WI 
o2o0000 


eQoooasac 
At OOM AO 
moves Or 
DU aad entry 
DNAS 
WOOaNSL 
eo ee @ @ @ 


ooo0000 


OTONTN 
DQoVYUOO 


Cocjooa 
xe—OnN Ot 
TMP NUON 
oT ON 
DST tet he 
ater INO) 
eeee% # 


ocoo00o 


(NM FuNO 
eg? Ne eee ee ee eee 


ROKK OK 


™ 
ph 


Dwr wees 


ooD0000 


oooceodo 
TOONS 
D aaah 0 
AONAOs 


MAN MAY Oat 
oo0000o 


Seer 


WOM 
WOW DO O 
NG at OF 
oe ee ee 
OCo0000 
{ Jo ug 


MW Tm 
Coo°oo° 


OoOoqooos 
Dean 
Ore Oey 
MOM™AOS 
DuUMWNOLN 
ODA Das Mtv 
ee ee ee ee 


oococec 


OTOININ 
YVUOQOoOvuUw 


eaoeaqge 
Oo? vNam 
Drm Wn 
NO Qe 
Ot ata 
NT Ue ant 
Ce ee ee 
ooo0eeco 

1 { 


NN FUN 


et Se? egg Se ee ee 


KO OK OOK 


ey 
~~ 


PUIG ON ANI CO ae 
oooo°oo 


eoeaede 
rom MS 
OH .rOmMmaoe 
SF RA wh 
WOMNDOO 
NAOH 
ee 8 © @ 
ooo00o 
I ots 


muNOorme 
oOo000c 


eEcssoe 
UVR ON OO ot 
MP OOrn 
TOW. Ory 
aN OO sy 
ON LD amt ery 
eee & © 8 


ooooce 


at OO wu OO 
fe? eee ee ee? eee ee 


me OK EO OK 


Tt 
ph 


ore seer 
Sooese 


oooocs 


NAOOAN 
am cam 


TAREQ CAA i et 
oo°o°eoco 


mK MK OK 


Ww 
fh 


Orr orm Ss 
coooeoec 


esaqgeoseG 
“FTN 
PWT © 
L$TOoOUomMmMsre 
Ota Oonm 
NOFNWIANM 
* ee © @¢ @ 


OoO000D 


FY ESCA CO ant 
oo0ooco 


SS Teas 
OP ont Fh OD 
KOoror 
mm Or Fo 
SNuWNOM 
a Mat MO FT OY 
eee 8e¢ @ @ 
ocooceo 
! titel 


OS ON Sa 
MWOUWVUOSO 


O2nD90005 
Oy med CS oh) ot CD 
CrONNEN 
NIN UN 


KK OK OK OS 


3 
ft 


Dormer sres 
ooo00°ec 


QOOQ0080 
MMH MOON 
DAnuranr 
NOWMAWNO 
Moon 
NOSFNOM,, 
°° © © @ @ 8 


oOO0O07D00 


DO oF NN ge on 
VOVUOU80 


Co00o0d0 
r~OrmmMoo 
Me moeor 
SOU Ch ae 
ante r~ Om 
WMO MT 
ere? 


oco0o°eo 
' \ 


at QUO PUNO 


me KOK OK OOK 


tht 
tf 


ores Ss 
oCoo0Co°O 


oooceaa 
hee Gantry 
NOt erow 
TNH OW 
mnt at OD ot 
NOIFNOM 
oe eo e @ © 


oCo000°0 


OYA CON ee 
ooooo0o 


Seoaceagae 
CP cme Pee free U7 
PNT OM © 
Oyoret 
NSP tet 
NO OO eet OO) AF od 
>. 2e¢e?e 6 6@ @ 
oCoocooo 
f toe 


MN Orh wn 
ooo00o0o 


eegege 
Phe OU amo 
DMVWAM Mo 
OM NID ant men 
a=at OO om 
ant OF UNAS CO 
eee 8 8 @ 


oOoog000 


OTN TN 
WwOWOVUO 


COoojsy00 
COOUOmeo 
UV WA 
OMmurn ON 
at FO a OC 
NIT Be OU 
® @© ee? @e@ @ 


oOooo00°0 
' t 


Pie ony > a, = 


NO FuUNO 


mC KOK OK OK OK 


@ 
t~ 


Owrh wh + 
ooo0o°0oo 


Oooooco 
DOnTruo 
TWODMOOD 
MODOMO 
NWO OOo 
NOITNONM, 
oe ee © 6 


oo00o00°o 


DOT WN 
Oooo0oco 


Qoo00o0o0ma 
TOA OM 
Over ory 
Oonwws ft 
mnt ot LAD OO) 
eevee 


oo00co 


OFMONTA 
Ovoveo 


oocoaa 
NNO 
MWe" FM 
NFIMOS 
NITNDOM 
NAW Pe OF 
ee © @ @ @ 


oOo00000 


in Om ey oh OT 


at OPIS 
ed 


<< KOK OOK 


oO 
t~ 


wo e* wee ae mt eet Bene Yo a eer ee, Oe eee ee 


70 





MEAN OF ERROR VAR CF ERROR 


VAR OF TRACK 


MEAN OF TRACK 


PAMe 


At IO NIM aw 
ooo0o0o0oo 


alefjeers 
MIN O Om © 
at OU ae 
CR OUP Pee cet ot 
DON — 
NPN NO 
oeeee 
Coooo0co 
| 


Mm OOS Nh 
oooo00°o 


aogaacae 
Waa m- 
aa OY SIA 
mayen Oo 
Not ODP oe 
Dats? OM 
eeee8t ¢ @ 


ooo 000 


OP MNS 
ovoove 


OoOoonoeo 
aTorrwy 
Us mt POUND 
ATO OW 
NITOOOM™ 
NM 2am 
eee ¢ @ @ 


Cooo°0e°o 


ey ey ey aly, oe 


KNOPF UIVO 


HE HE HK OK OK OK 


o 
© 


- 


Ser rer 
oooo7to 


oOoooao0 
DON om 
DROS 
PN at et oN 
eI OOSCWN 
Manon 
© ¢@¢ @ @ 8 


oooooo 


NNOWNM = 
Cooecos 


CcCocoom 


PUN OLIN in 
OoCoo°co 


ogee oes 
aor~on™ 
gowermDuon 
mt Orn 
mNOOM 
COD ST et me 
* ¢ @ 0@ @ @ 


Oooo00o 


OS MN Sm 
Ovve09o 


OOOocoo 
at Det AUS 
NROODeN 
OD ON et et 
ONS Poet LN 
QUIS OO OO 
ee ¢ ¢ @ © ¢@ 
200000 

' ' 


(ae fn fan gi, A 


AN TUN 


MK MOK MOK 


~_ 
eo 


WN SPUN 
OoOoo°ooc 


NIN at OY ON oe 
oooeco°o 


oo qeaas 
HONNOr 
PN OUNO AN 
ODD Oty 
rh oOMm 
AAT 

oe @ © @ @ 
oooooo 

‘4 


mA OU 
ooo0o0°0o 


Seooaes 
™ MOE-T DOO 
Om OOM 
WN D> AU ow 
DBOonmMownr 
PLAT et sal 
eee ee 


oooo7do 


OTMAM SOY 
VOOCOVLY 


monrnona 
TOOUONEe 
O'Nnt D-F 
NOS MO 
$A) Ory Pe mt OO 
ONY TY ET ant OM ee 
* ©@ @ © &¢ 


OOoonDco 
I ' 


om a om oe 


MANTIS UNO 


eet Se eee Ce Se 


HE HE HE OC OOK 


N 
@ 


LN AS UNT M 
Oooqoo7Te 


ooooocr 
KOON O 
OOO Oe 
Or @On~wn 
Mm MLN co me 
Om NOUN 
ee. ee¢ ¢ 


OoDdo0oog 


NN et ON 
oooo°o°o 


SeooqgeoS 
Paw Th OT 
KON S 
SNM HN 
COON Oem hm wr 
gan (tN mE TN) omet 
e¢ ee © @ ¢@ 


ooo0o°o0°o 


MmINOIA™ WL 
oooo0o°o 


ooooeS 
mom A 
INO NICOWY 
WO Oat eym 
Pet SHON 
aA OO 
ose ee 


ooooec 


Pty ly en ee ly 


a NTT UNO 


wr tee Cee SOP te ee 


MME HE OE OHO 


mM 
eo 


WPUIAT MH 
ohaloloal ole) 


Soooo0oa 
att Fat OT 
MUNA OS 
Or OO 
OFT MOKWM 
MAO TNIN 
¢ ¢ e@ © © # 


eo leloleleale) 


NIN NO ION 
ooo0ooeeo 


Saqgjges 
NILA OF & oO 
NO™oOmMO 
BONN 
aANOIOOC 
Ne RO oP ee 
ee © © @ 8 


Oooo0coo 
\ 


DUN OLIN 
oooo0ceo 


o90000G 
=2O0™r 
PIAA SV 
LONHDS 
rT Omit 
Cm mer un 
ee e8 @ @ 


Oooooo 


OTMOM TA 
VWOoOOVOCOS 


OQOOuo0o 
uIsuIONeas 
DBYVNO eS 
~~ OOM 
MAO SUNN 
ONS ON nt OH et 
eee 6 


oooocofo 
' { 


Se ey ee ee ee 
Haruno 


HC HK HE OK HOOK 


wT 
wo 


UNS LF Pe wy 
oooooo 


Oooeonao 
NNN O = 
OD OOq 
OH Fh AO 
ONDA mr 
ONT ton 

o@eef © @ 


Ooo000 


OFNMAMSIO 
OVCOVOS 


COD0D8 
ah oN 
wor vO 
NOAA om 
MAN TOME 
CIO Mt OT 
os 6 © @ 


oooo0o 
J 


amt TLD 


Rene Cane? tee tee Se See 


2S HK HK OK OKO 


uw 
@ 


a era ance a a TO US APN AMRIT S ee mm rm St ont 


a eee 


ae 


TE Se EY SO Uae TENS Go SF ET 


Ty 





AE PENUTA Te 


Meee olSs OF ATIGS SIMULATION WITH POSITION RESET AND KALMAN 
FPILTERING 


72 









’ 3s 
_ oe a3 aa 










Ot oe 
2296 ODFTIEOS ATTR iat tT 0th 40 | 
~——  grtserzT - : 





Hd 


Se ESOS OO ey eA are at ot Se ee el SO OOS eh SN ee Oct AEs i se 





UH? a dh did Ald TE (S8€) 


Figure 21 - SQRT OF DOWN RANGE VARIANCE CATIGS WITH 
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Figure 23 - SQRTe OF DOWN RANGE VELOCITY VARIANCE CATIGS 
WITH POSITION RESET.AND KALMAN FILTERING? 


wo 





FI/SEC 


1d Ah A 


— 


——— 


bi 


4d 


See 


4d 


i 


ee 


an | 


Ha Hi? itis iti Hf AID TIME (EC) 


Figure 24 - SQPTe OF CROSS-RANGE VELOCITY VARIANCE CATIGS 
WITH POSTION RESET AND KALMAN FILTERING? 
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Figure 25 - SQRTe OF DOWN-RANGE VARIANCE CATIGS WITH POSIT 
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78 





S 


= = 


= eT A eS ES 


Hil 2 a ih dT Ald TIME (SEC) 


mequiee 26 —= SORTS OF CROSS-RANGS VARIANCE SATIGS WITH POSIT 
RESET AND KALMAN FILTERING AT THE X5 NOISE LEVEL? 


We 





i 


: 
TH 


‘dE 


<a 


Wt 
a0 


Wi 


ni U2 a (ith Uf WIA TYRE (560) 


oa -—— =e — 2. 8 ee te ee ee tr ee me ee eget 


mies —- SORTS OF DOWN-RANGE VELOCITY VARIANCE CATIGS 
Peet POlog) RESET AND KALMAN FILTERING AT THE x5 NOISE LEVEL? 


80 





. ce ae ne ee eg 


FT/SEC 
Md 
2 
Hd 
A] 
52 
H J 


2H 


Figure 28 - 


i Af ail Aid TAME (Ot) 


SQRTe OF CROSS-RANGE VELOCITY VARIANCE CATIGS 


WITH POSIT RESET AND KALMAN FILTERING AT THE X5 NOISE LEVEL? 


8 1 





PARTIAL LISTING 


X (i,j) 
THETA (i,j) 
Wor (1, >) 
META (1,3) 
MEER (i, 4) 
XBVAR (i,3) 
XMEAN (i, 3) 
XVAR (i, 3) 
A (J) 

AWXY 
AWYY 
AWYZ 
AWZZ 
AWXX 


DEL VX 
DELVY 


APPENDIX E 


OF SYMBOLS AND NOMENCLATURE OF SIMULATION 


PROGRAM 


i-th missile state at time j 


i-th angular state variable at time j 


estimated i-th state at time j 


estimated i-th angular state at time j 


mean of i-th state at time j 


variance of i-th state at time j 


error mean of i-th estimated state 


error variance of i-th estimated state 


thrust acceleration at time j 


velocity changes due to angular rotation 


changes in wind components 
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BETA1 
BETA2 
BETA3 
YLIX 
WiyY 
YIZ 


XPOS 
n2OS 


G1 
G2 


OM EGAX 
OMEGAY 


XINT 


SL 


XFIN (i) 


NTERM 


VXO 


VYO 


Ix 


DITAHTAX 


dthtxn 


Z2NI (3) 


body referenced accelerations 


inertial referenced accelerations 


micrad sensed positions 


Kalman gains for filter 


extra states for Kalman filter 


dummy variable for output 


Change in drift angle 


final vaiue of i-th state per track 


maximum time steps allowed 


velocity of wind down-range 


velocity of wind cross-range 


seed number tor random number generators 


change in thetax 


measured change in thetax 


total tracks through time j 
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XIFIN (i) 
track 


XBFIN (i) 


XIBEN (i) 


XBF V (2) 


NNA 


ITENSB 


SIG EO 


SIGW 


SaG K 


SIGEG 


SIGKG 
factor 


SIGT 
theta 


final value of estimated i-th state per 


mean of XFIN(i) 


mean of XIFIN(i) 


variance of XFIN(i) 


humber of gyros Simulated 


number of accelerometers Simulated 


size of ensemble 


std.deviation 


std.deviation 


std.deviation 


std.deviation 


of 


of 


of 


std.deviation 


std.deviation 
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gyro bias 


gyro random walk 


gyro scale factor 


accelerometer bias 


of accelerometer scale 
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CIV? 


AAAQ 


AQ 


AAAA AAQ AAA 


i 


18 


Ve 


20 


30 


50 


> 1 


GENERATE GYRO BIAS "EO" 
CALL SNORM (IX,EO,N) 

DO 18 KN=1,N 

EC (KN) =SIGEO* £0 (KN) 


GENERATE THE RANDOM SCALE FACTOR "Zk" 
CALL SNORM (EX, AIS 
DO 19 KN= 


ZK (KN) = S SIG KEZK (KN) 
GENERATE RANDOM INPUTS TO ACCELEROMETER PACKAGE 


GENERATE ACCELEROMETER BIAS "FOG" 
CALL SNORM (ZX (BOG, NWA) 

DO 20 KN=1,NNA 
BOG (KN) =SIGEG*EOG (KN) 


GENERATE ACCELEROMETER SCALE FACTOR "KG" 
(IX,WG, NNA) 

NNA 

reas 


N 
eee 


HAS 

CW O 

mz WH Il bg 

ERG GQ HO 


ONO ES COR) err "KC CNY 
~y"“es es 8 83 Ht Uo lme ton xt 


~ 
oO 


( 
( 
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Pstg 
Cy 
mG ort res 


Oto th tee EEE te tr ON OM 
emer naerne @ SM thm ONY OV iets OUmCUN 
HUN WU dao Grree s £s 


NONE WN etd tit llLNN KRDO OW Ee eer 
ro PS es es es PS OOO rm = Gir 


a a rT a a or) 
Do bg Do ed Bs ed FESS Cy EY De es 


++ e+ 4+ 406) 


M IS THE NOISE INPUT FOR EACH ACCELEROMETER 
1) =E0G (1 ep ) xa (J) 

2) =EOG (2) +WG (2 

3) =E0G (3 

ERATE THE BIAS TERM DUE TO RANDOM WALK 
SNORM (IX,G,N) 

0 JI=1,N 

) =SIGW*G (JI) 

I 


S THE CHANGE IN THE ANGLE BETWEEN THE COMPUT 
RDINATE PLANS AND THE ACTUAL COORDINATE PLANE 


DO 51 Ji=1, 

BSI (JI) =£0 (J 5I) #6 (37) 
DELVX=VXO-VX 
DELVY=VYO-VY 


SOmruUne INERTIAL POSITION 


Bee EOS Dae 1) CORNER YEE (2 
DTHTYM=DTHTAY+PSTI (2) +DTHTAY*ZK (2 


QOQ21GQ) QQ Peo P< Od OG RIES De Pe te Be PH et cI OHOOUNO aH 


OH G 


Q@) 
= (2 .f& 


re 
ROE 


BD 


QWe |OaONM 
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AAA AAAQN 


QAAQ 


AAANA 


AAAAQANQAAQNAAQ AAQ AADQ AAAAYA 


QAAX1AND 


DTHTZM=DTHTAZ+PSI (3) +DTHTAZ*ZK (3) 


THE GAINS G1 AND G2 ARE THE KALMAN GAINS GENERATED AS 
EUNGTION OF TIME 


ro 503) 
(Oe (tra X+1 


COMMANDED HEADING CHANGE IS SUBTRACTED FROM THE 
ERVED HEADING CHANGE 


J=DTHTYM-DTHTAY 
J=DTHTXM-DTHTAX 


EQUATIONS FOLLOW 

1,9) #G1* (DLXI- OMEGAX) 
LXJ-OMEG 

385) GE (D Y J-OMEGAY) 

LY J-OMEGAY) 

S ARE USED TO PREDICT THE NEXT 
TOR 

i {3,3} SOEMEAYSOneGhY 

I 


=". Ts 


E 
S 
- 
x 
iE 
E 
E 
E 
E 
E 
A 


2,J) tDTHTAY+OMEGAY 
So) Fo LET 


aoe =o OnOd = 


rgn a & 


PHO i a 
Sodd me gar 


TIMEX=1 
SENSED ACCELEROMETER INPUTS IN THE BODY AXIS FRAME 
IN THE X(BODY FRAME) DIRECTION 

BETAI=A (J) -DELVX+DELVY*TYETA (2,J) +GAM (1) 

IN THE Y(BODY FRAME) DIRECTION 
BETA2=DELVX*THETA (2,d) +DELVY+GAM (2) 

IN THE VERTICAL(BODY FRAME) DIRECTION 


BETA3=0.0 

THE PURPOSE OF THIS SECTION IS TO GENERATE THE 
INERTIAL ESTIMATES OF POSITICN BASED ON A PURE 
INERTIAL COMPUTATION 

AWXXI IS THE ACCELERATIONS DUE TO HEADING CHANGE 
AFFECTING THE X DIRECTION FROM THE ANGLE CHANGES THETA 
X. SIMILARLY AWXYL IS THE ACCELERATIONS AFFECTING 
THE X DIRECTION DUE TO THETAY 

AWXXT=XI (6 J) *DTATAX 

AWKXYI=-X 4 J) *DTHTAY 

AWYYI=X1(2,3) *DTHTAY 

AWYZI=X1L (6,0) *DTHTZM 

AWZZiI=- X hyd) SDRETZE 

AWZXI=-X1(2,J) *DTHTAX 

DTHTXM=0.0 

DTHTYM=0.0 

DTHTZM=0.0 

DTHTAX=0.0 

DTHTAY=0.90 

DTHTAZ=0.0 

SENSED ACCELERATIONS IN THE BODY FRAME ARE CONVERTED 
TO THE INERTIAL FRAME. 

YIX=BETA1-BETA2*THETAL (2,3) 
YIY=BETA1*THETAL (2,0) +BETA2 

IF (J.FQ.15)GO TO 48 

IF (J«NE.80)GO TO 49 
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GENERATE THE RANDOM ERROR IN THE POSITION MEASUREMENT 
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